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Combustion  Zone  Characterization  of  G02/GH2  Rocket  Using  Laser-Induced 

Fluorescence  of  OH. 
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Propulsion  Engineering  Research  Center 
and 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

With  recent  interest  in  gas/gas  injectors  for  use  in  rocket  combustors,  there  is  a 
critical  need  for  experiments  that  address  this  combustion  process  in  terms  of  detailed 
flowfield  measurements.  Such  measurements  would  also  serve  as  a data  base  for 
validating  computational  fluid  dynamic  (CFD)  computer  codes.  A series  of  studies  have 
been  undertaken  at  the  Propulsion  Engineering  Research  Center  (PERC)  at  the 
Pennsylvania  State  University  to  measure  various  parameters  such  as  velocity,  species 
concentration , and  temperature  downstream  of  a shear  coaxial  injector  in  an  optically 
accessible  uni-element  rocket  chamber.  Techniques  applied  to  this  study  to  date  include 
the  following:  laser  Doppler  velocimetry  (LDV)  for  velocity;  laser  light  scattering  (LLS) 
for  flow  visualization  and  estimating  mixture  fraction  and  density;  laser-induced 
fluorescence  (LIF)  of  hydroxyl  radicals  (OH)  to  determine  the  characteristics  and  extent 
of  the  reaction  zone;  and  Raman  spectroscopy  to  measure  major  species  concentrations 
and  temperature.  The  results  of  the  LIF  studies  are  presented  here. 

The  OH  molecule  is  a key  intermediate  in  hydrocarbon  and  hydrogen  combustion. 
High  OH  concentration,  indicated  by  high  fluorescence  intensity,  mark  the  location  of  the 
primary  reaction  zone  where  the  oxidizer  to  fuel  ratio  is  nearly  stoichiometric.  Two- 
dimensional  imaging  of  LIF  near  the  injector  face  provides  a qualitative  view  of  the 
reaction  zone  structure.  Two-dimensional  LIF  was  limited  to  qualitative  measurements 
near  the  injector  face  due  to  poor  signal  to  noise  ratio  with  the  present  experimental  setup. 

One-dimensional  measurements  of  LIF,  which  provide  a radial  profile  of  relative 
OH  concentration,  have  been  made  at  several  axial  locations  in  the  combustion  chamber. 
Results  from  multiple  images,  typically  120,  have  been  averaged  to  yield  average  OH 
profiles  at  each  axial  location  probed.  Probability  density  functions  (PDF)  of  OH  peak 
widths  and  locations  show  that  the  reaction  zone  is  thin  near  the  injector  face  as  expected 
and  remains  thin  as  the  flow  progresses  downstream.  Also,  the  increase  in  widths  of  the 
average  OH  peaks  as  the  flow  progresses  downstream  is  due  to  movement  of  the  thin 
reaction  zone  rather  than  an  increase  in  individual  OH  peaks.  This  analysis  indicates  the 
flame  is  a wrinkled  laminar  flame  front  in  the  region  probed. 
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■Iff  MOTIVATION 


Need  For  Low  Cost,  Reusable,  and  Reliable 
Propulsion  Systems  Require  Detailed 
Understanding  Of  Combustion  Phenomena 

• Full  Scale  Tests  Are  Expensive 

• Computer  Technology  Enables  Detailed  Modeling 

• Advances  In  Optical  Diagnostic  Techniques 
Enhance  Measurement  Capabilities 
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Obtain  Data  Where  Boundary  Conditions  Are  Well 
Specified  And  Provide  To  Rocket  Research 
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Indicates  Reaction  Zone  Location  and  Structure 


■IV  LASER-INDUCED 

FLUORESCENCE  OF  OH 
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CED 

FLUORESCENCE 

Experimental  Setup 
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■ B PEAK  POSITION  PDFs 
LIF  PROFILES 


■Ilf  SUMMARY  A CONCLUSIONS 
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Average  LIF  Profiles  Useful  for  General  Profile  Shape  and 
Location  Comparisons 

Velocity  and  Light  Scattering  Data  Also  Available 


Spatially  Resolved  Species  Measurements  in  a G02/GH2  Propellant  Rocket  £ 5 

M.J.  Foust,  T.  Ni  and  R.J.  Santoro 

Propulsion  Engineering  Research  Center 
and 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Statement  of  Problem 

Local  species  concentration  and  temperature  are  among  the  most  important  parameters  to 
characterize  in  a combustion  system.  However,  the  harsh  environment  in  rocket  chambers  under  hot-fire 
conditions  limits  the  applicability  of  conventional  probing  methods  for  the  acquisition  of  this  information. 
Laser-based  diagnostic  methods  show  great  promise  for  achieving  this  goal  by  providing  instantaneous 
images  of  species  concentrations  and  temperature. 

Objective  of  Work 

The  objective  of  the  current  work  is  to  develop  a non-intrusive  technique  to  experimentally 
determine  the  major  species  and  temperature  field  in  the  combustion  chamber  of  an  uni-element  rocket 
for  a G02/GH2  propellant  combination. 

Approach 

The  experiments  were  conducted  at  the  Cryogenic  Combustion  Laboratory  at  Penn  State 
University.  The  uni-element  rocket  chamber  used  is  modular  in  design  and  is  easily  configured  to 
provide  optical  access  along  the  chamber  length.  A shear  coaxial  injector  was  used  to  introduce  G02  and 
GH2  into  the  combustion  chamber.  The  nominal  flow  rates  of  G02  and  GH2  are  0. 1 Ib/s  and  0.025  Ib/s, 
respectively,  resulting  in  an  O/F  ratio  of  four.  The  experiments  were  for  a chamber  pressure  of  190  psia. 

One-dimensional  profiles  of  species  concentrations  and  temperature  were  measured  by  using  laser- 
induced  spontaneous  Raman  spectroscopy.  The  Raman  system  consists  of  a flash  pumped  dye  laser 
operating  at  10  Hz  and  an  intensified  CCD  camera.  The  dye  laser  has  a typical  pulse  energy  of  2 J at 
511  nm  and  a pulse  duration  of  5 /xs.  The  Raman  emission  was  detected  at  a right  angle  to  the  laser 
beam.  A narrow  band  interference  filter  was  placed  in  front  of  the  camera  to  selectively  measure  the 
number  density  of  the  species.  For  each  rocket  firing,  50  single-shot  Raman  images  and  50  background 
images  were  captured.  By  using  different  optical  filters,  Raman  images  of  oxygen,  hydrogen,  water  and 
nitrogen  (used  for  window  cooling)  were  obtained.  Measurements  were  conducted  at  1,  2,  and  5 inch 
downstream  of  the  injector  face.  The  ratio  of  the  signal  to  background  level  for  hydrogen  and  oxygen 
Raman  images  at  1 inch  downstream  is  about  10.  Further  downstream,  the  background  luminosity 
increases  significantly.  Thus,  the  species  concentration  can  only  be  determined  from  averaged  Raman 
images.  The  temperature  profiles  were  calculated  from  averaged  data  of  total  species  number  density 
using  the  ideal  gas  law.  Since  the  Raman  signal  is  stronger  in  lower  temperature  regions,  the  averaged 
temperature  generally  underestimates  the  temperature  in  regions  where  temperature  fluctuates  highly. 

Conclusions 

Single-shot  and  averaged  profiles  of  species  concentration  have  been  measured  under  combusting 
conditions.  These  results  demonstrate  that  the  laser-based  technique  can  be  effectively  applied  for  in-situ 
measurements  in  a rocket  chamber.  Experiments  with  an  improved  detection  system  for  obtaining  images 
of  instantaneous  and  simultaneous  multi-species  concentration  and  temperature  are  underway. 
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go2/gh2  rocket  test  conditions 

Shear  Coaxial  Injector 
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EXPERIMENTAL  SETUP 


• Optimized  System  for  Collecting  Weak  Raman  Signal 
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EXPERIMENTAL  MEASUREMENTS 

Mean  Density  Comparisons  @ 25.4  mm  from  Injector  f ace 
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Oxygen  Mole  Fraction  Comparisons 
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SIMULTANEOUS  SPECIES  MEASUREMENTS 

25 A mm  from  Injector  Face 

Log  Intensity  Scale  Linear  Intensity  Scale 
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Fluctuating  Pressure  Analysis  of  a 2-D  SSME  Nozzle  Air  Flow  Test 


Darren  Reed 

Induced  Environments  Branch 

NASA  Marshall  Space  Flight  Center,  Huntsville,  Alabama  35812 


■ 


c;y  a. 


Homero  Hidalgo 
SSME  Assurance  Office 

NAS  A Marshall  Space  Flight  Center,  Huntsville,  Alabama  35812 

Abstract 

To  better  understand  the  SSME  startup/shutdown  transients,  an  airflow  test  of  a 2-D  nozzle  was 
conducted  at  MSFC’s  Trisonic  wind  tunnel.  Photographic  and  other  instrumentation  show 
during  a SSME  start  large  nozzle  shell  distortions  occur  as  the  Mach  disk  is  passing  through  the 
nozzle.  During  the  earlier  development  of  the  SSME,  this  startup  transient  resulted  in  a low 
cycle  fatigue  failure  of  one  of  the  LH2  feedlines.  The  2-D  SSME  nozzle  test  was  designed  to 
measure  the  static  and  fluctuating  pressure  environment  and  color  schlieren  video  during  the 
startup  and  shutdown  phases  of  the  run  profile. 

The  model  consisted  of  two  identical  blocks  having  the  same  inner  contour  of  the  SSME  nozzle. 
The  sides  of  the  nozzle  were  made  of  glass  for  schlieren  photography.  The  upper  block  was 
instrumented  for  static  pressure  measurements.  The  lower  block  was  instrumented  with  thirteen 
Entron  fluctuating  pressure  transducers.  Steady  state  and  slow  sweep  flows  were  tested  for  three 
back  pressure  conditions  (0.5-2.0  psi,  7 psi,  14  psi.)  The  static  pressure  data  was  acquired  by  a 
scanning  pressure  system.  The  fluctuating  pressure  data  was  recorded  onto  a VHS  analog  tape 
recorder.  The  video,  static  pressure,  and  fluctuating  pressure  data  were  time  synchronized  for 
data  correlation. 


The  shlieren  video  clearly  shows  a lambda  (k)  shock  foot  moving  down  the  throat  during  the 
slow  sweep.  The  fluctuating  pressure  RMS  time  histories  show  the  levels  increase  as  the 
downstream  foot  of  the  lambda  shock  approaches.  When  the  shock  foot  is  directly  above  the 
transducer,  levels  decrease  about  50%.  When  the  upstream  leg  of  the  lambda  shock  approaches 
the  transducer  the  level  quickly  jumps  up  to  twice  the  downstream  leg  values.  After  the  upstream 
leg  of  the  lambda  shock  passes  the  transducer,  the  level  falls  down  to  the  noise  floor  of  the 
measurement. 

Schlieren  video,  model  configuration,  fluctuating  pressure  time  histories,  power  spectrum 
densities  of  the  test  will  be  shown.  Future  2-D  nozzle  tests  and  plans  for  a 3-D  nozzle  facility 
will  be  addressed. 
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incorrect  weld  material 

resolved  by  adding  nickel  plating  to  tee  weld  joints,  added  steam  loop  to  coolant  line 


726 


)r.  in  CD 
® in  co 
c cr  a 
o o uj 

*r  o m 
■C  O)  Q) 

3 a ® 
rr  -a  tr 


a>  x g 
8 2 t 
a.  ® g 
to  E u 
_ o 
n I 
£ 

CO 

L. 

(0 


Q 

■ 

CM 

CO 

«4—  4-J 

o w 

!•§ 
so  5 
C LL 

< 


□ 

(/> 

W 

® 


< 

0) 

N 

N 

o 


£L  Z 

m UJ 

c § 
’■5  co 
§ co 


m 

© 

> 

o 

© 

8 BBBIWIB^ 

JD 

o 

% 

© 

h“ 


0 

E c 
5 2 

Q.  CO 
X CO 

® ^ 

0 0 

c N 
C N 
3 O 
C 

■5  Hi 

iS 

co  CO 

CO  S 

5 «- 

5=  O 

5 0 
N "O 
N O 

■o'  o 
re  a 
0 c 
o 


o 


re 

o 


0 
■*— < 
4-> 

0 

13 

O 


>*  - O 

ss  c 

.ram  CC 

o „ 

5 8 

0 
c 
c 

3 
h- 

■o 

c 


co 

c 

3 

0 g 

JC 

T3  Q. 
C w 
0 Q 

•4-  8 

£ 

0 -a 
13  0 
c 
3 


o 

o 


£0 

0 

O) 

c 

‘to 

3 


CO 

0 


3 

CO 

CO 

0 

Ira. 

Q. 

03 


o 

’E 

o 

CO 


0 

3 

*■* 

o 

3 


_ — ^ 
H O 


O 

c 


0 
0 

* £ 
^ 3 
0 «o 
0 

O 0 


LL 

CO 


0 

X3 

0 
.*— « 
O 
3 
"O 


_ 0 

f 2L 

2 <s 

■g  § 

= •! 

0 > 

is 


+5  0 
0 = 
c x: 
— o 

©00 

Stj 

> 0 


0 
£ 

0 0 3 

O O 


TO 

C 

*3 

3 

■a 

0 

N 

N 

o 

c 

0 

£ 0 
Z g 
o > 

■5  0 
of 

§® 
o c 

Era 
- E 
0 0 

CO  -C 

3 ° 

o? 


0 « 

^ i 

0 5 

o c 

W o 

© "O 
JZ  ^ 

4-<  3 

•«-  x: 
o 0 

© 05 
0 C 
■O  -3 

11 
© ^ 
= O 
X 0 
O A 

?S 

■a  a 

o5 

0 

0 <S 
QC  0 


■a 

c 

ra 

0 

0 0 
^ > 
a 0 
o — ■ 
jc  >* 
0 TJ 


0 


ra 
0 

H—  C 

O 3 

£ © 
o>£ 

0 0 

£ 3 

0 0 

© ra 

■■§  1 

1 2 

« 3 

£ 0 

3 

0 

c g> 

M=  0 

0 > 

■O  5 

a £ 

1 § 

3 
■o 

8 0 

ra  g- 
* ro 


w 

X 

33 

3 

o 


i 2 

4-1 


0 


0 0 


© 3 E 
Q.  CO  3 

0 0 J= 

i=  2 a 
3 CL  0 

0 z:  o. 
0 0)0 
0 c 
■—  0 

Q-  ro£ 
0 o t! 

+*  o > 


S 3 


0 0: 
0 0 


— x: 


0 

0 


SZ  SZ 

h*  I— 


© O 
O 


-*V  A 


727 


•Sk  in  co 

3 jp  % 

£ ll  □ 
®QLU 

o ^ 

t o "O 
£ O)  Q) 
.2»«  $ 
E S c£ 

« x s 

o o £ 
Q.  5 to 
w £ Q 

= O 

*5  ^ 

jc 

m 

Xm 

m 


o m 

« »2 

£g 

C0  ii 
C LL 
< 

0 


< 

0) 

N 

N 

O 


m 
m 

0 
l. 

Q.  Z 

m ^ 

c SE 

1 m 
§ CO 

o 

3 

nJLi 


CD 

C 

O 

Q, 

□ : z: 

o 

m 

m 

a 

% 

0 

H 

13 

C 

CO 

0 

"D 

O 


0 

Q. 

CQ 

W 

SU 

3 

O 


o 

0 

HI 
2 
m 
m 

1 

cm 


x 

© 


(0 

o 


(0 

© 

JC 

o 

c 

03 

C£^ 

© 

© 


W 

© 


CO 

00 

II 

g 

"■«-> 

CO 

CC 

© 

© 


•==  £ 

- q 
II  w 

B 11 

05  B 
C T3 

© > 
-J  > 

© © 

N N 
N N 
O O 


© 

© 

sz 

o 

c 

0 
in 
ii 

4-s 

£ 

05 

© 

1 

44 

'x 

LU 

© 

N 

N 

O 


© 

© 

O 

c 

00 

to 

in 

■ 

o 

ii 


05 

’5 

x 


© 

o 


JC 

o 

o 

CO 

i- 

© 

5! 

o 


■O 

© 

■O 

o 

o 

© 

X 

0) 

© 

o 

3 

■o 

CO 

c 

© 

H 

© 

3 

CO 

CO 

© 


JC 

o 

o 

CQ 

a, 

o 


O)  CO 

.£  s 

15  o 

I 0. 

o 

3 
LL 
CM 


■—  4-4  m 

♦J  O L 

Co 

C 
0 

E 

3 


CO 

c 


I 

© 

r 

o 

0. 


Z.  t=  r '-5  2 


3 

© 

W 

© 

v. 

X 

o 

*3 

© 

4-> 

0 


X 

su 

© 

N 

N 

O 


© 

© 

1_ 

3 

. W 
© 
© 

CL 

O 

4-J 

© 

4-» 

0 

TJ 

C 

© 

©" 


© 

© 

Q. 

E 

© 

l- 

cd  « 


E 

0 

E 


© 

*-» 

o 


r O 


5 2® 


_ 3 2 

3 © > 

CD  W 


s £ i 


< Z X 2 |= 


^ A A A 


■C  — oo  00 


0 

■o 

o 


■ © 

^ © £ 

S’  o o 

■—  H-  0 

°0  _ 
€0  ~ ~ 


0} 

c 

o 

«n» 

4*5 

omn 

TJ 

C 

o 

o 


o 

© 

JC 

CL 

OT 

O 

£ 

4-i 

© 

T2 

C 

© 

t© 

'© 

Q. 

h» 

© 

© 

Q. 

CM 

« 

C 

© 

44 

‘“5 

c 

o 

o 

© 

1- 

3 

OT 

OT 

© 


X 

© 

© 

N 

N 

O 

z 

CO 


OT  A „ 

fl)  * * * 

h" 


728 


Slow  sweep  runs 

5 steady  state  runs  at  predetermined  shock  locations 
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EXPERIMENTAL  INVESTIGATION  OF  THE  EFFECTS  OF  ACCELERATION  ON 
HEAT  TRANSFER  IN  THE  TURBULENT  BOUNDARY  LAYER 


Walid  M.  Chakroun 
and 

Robert  P.  Taylor 


Department  of  Mechanical  Engineering 
Mississippi  State  University 
Mississippi  State,  MS  39762 


Abstract 

Rocket  propulsion  system  components  such  as  turbines/pumps  and  nozzles  often  have 
aerodynamically-rough  surfaces  or  surfaces  which  become  rough  during  operation.  Also,  these 
surfaces  are  often  in  regions  of  accelerating  flow.  The  interaction  between  surface  roughness  and 
acceleration  is  complicated  and  not  predicted  by  a simple  superposition  of  flat-plate  rough-wall 
correlations  and  smooth- wall  acceleration  effects. 

For  the  smooth  wall,  acceleration  causes  a decrease  in  the  Stanton  number  when  compared 
with  equivalent  unaccelerated  flow.  When  the  acceleration  is  strong  enough,  the  turbulent 
boundary-layer  heat  transfer  rates  will  approach  those  of  a laminar  flow  and  the  boundary  layer 
is  said  to  have  relaminarized.  Under  proper  conditions,  rough-wall  accelerated  flow  can  have  the 
opposite  behavior  with  increasing  Stanton  numbers  and  hence  much  larger  heat  transfer  rates. 

The  objective  of  this  research  was  to  experimentally  investigate  the  combined  effects  of 
freestream  acceleration  and  surface  roughness  on  heat  transfer  and  fluid  flow  in  the  turbulent 
boundary  layer.  The  experiments  included  a variety  of  flow  conditions  ranging  from 
aerodynamically-smooth  to  transitionally-rough  to  fully-rough  boundary  layers  with  accelerations 
ranging  from  moderate  to  moderately  strong.  The  test  surfaces  used  were  a smooth-wall  test 
surface  and  two  rough-wall  test  surfaces  which  were  roughened  with  1.27  mm  diameter 
hemispheres  spaced  2 and  4 base  diameters  apart  in  a staggered  array.  The  experiments  were 
conducted  in  the  Turbulent  Heat  Transfer  Test  Facility  in  the  mechanical  engineering  laboratories 
at  Mississippi  State  University.  The  measurements  consisted  of  Stanton  number  distributions, 
mean-temperature  profiles,  skin-friction  distributions,  mean-velocity  profiles,  turbulence-intensity 
profiles,  and  Reynolds-stress  profiles. 

The  Stanton  numbers  for  the  rough-wall  experiments  increased  with  acceleration.  For 
aerodynamically-smooth  and  transitionally-rough  boundary  layers,  the  effect  of  the  roughness  is 
not  seen  immediately  at  the  beginning  of  the  accelerated  region  as  it  is  for  fully-rough  boundary 
layers;  however,  as  the  boundary  layer  thins  under  acceleration,  the  surface  becomes  relatively 
rougher  resulting  in  a sharp  increase  in  Stanton  number. 
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numbers,  from  Coleman  (1976) 


In  a smooth-wall  boundary  layer,  acceleration 
stretches  the  eddies  reducing  the  trans- 
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OBJECTIVES  OF  THIS  WORK 
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Velocity  and  turbulence  profiles  measured  with  hot-wire 
anemometry. 
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Computational  and  Experimental  Efforts  € 

in  Gravity  Probe  B Microthruster  Analysis 

/ O Q/  f 

Gravity  Probe  B,  an  experiment  to  test  the  theory  of  relativity,  will  be 
launched  near  the  turn  of  the  millennium.  Due  to  the  precise  pointing 
requirements  needed  to  successfully  carry  out  this  experiment,  the  satellite 
will  use  sixteen  proportionally  controlled  micro  thrusters  as  a main 
component  of  the  attitude  control  system.  These  microthrusters  use  the 
helium  boil-off  from  the  on-board  dewar  as  propellant. 

Marshall  Space  Flight  Center,  overseeing  the  project,  verified  the  design  of 
the  thruster  flow  path  by  both  computational  and  experimental  methods. 

The  flow  performance  of  the  thruster  has  been  adequately  characterized. 

Graphs  show  specific  impulse,  thrust  coefficient,  discharge  coefficient,  and 
mass  flow  rate  trends.  Value  was  added  to  the  program  through  gained 
confidence  in  the  design  of  the  thruster  and  through  evaluation  of  some 
design  trade-offs. 

This  work  may  be  valuable  in  the  future  due  to  the  possible  need  of  small 
thrusters  on  spacecraft  that  have  precise  pointing  requirements. 
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Computational  and  Experimental  Efforts 
in  GPB  Microthruster  Analysis 


Computational  Fluid  Dynamics  Branch 
Fluid  Dynamics  Division 
Structures  and  Dynamics  Laboratory 
George  C.  Marshall  Space  Flight  Center 


Background 


• Gravity  Probe  B (GPB) 

- A satellite  borne  relativity  experiment 

- Requires  precise  pointing  control  and  acceleration  free 
environment  to  be  provided  by  attitude  control  system 

• Micro  thrusters 

- Helium  gas  from  dewar  boiloff  used  as  propellant 

- Sixteen  microthrusters  on  spacecraft;  used  for  orbit  trim,  spin-up, 
spin-down,  and  attitude  control 

- Concerns  about  mission  lifetime  and  control  saturation 
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use  of  large  liquid  helium  dewars 
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Method 


Direct  Simulation  Monte  Carlo 
(DSMC): 


Limits  on  areas  of  application 

Slow  - not  useful  for  large 
parametric  studies 

Works  well  for  low  Reynolds 
number  flows,  costly  to  use 
for  higher  Reynolds  number 
flows 

Gives  good  characterization  of  the 
flowfield 


Experiment: 


Covers  all  of  the  flowfield 

Very  fast  once  hardware  is  in  place 
useful  for  parametric  studies, 
useful  for  assessment  of 
configuration  change 

Possible  data  scatter  at  low 
Reynolds  numbers,  but 
works  well  for  higher 
Reynolds  number  flows 
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Results 

• Requirements: 


- Thrust  > 8 mN  at  Pinlet>  9.7torr  and  mass  flow  < 1.52E-05  lbm/s 

- Thrust  < 0.05  mN  at  Pinlet  < 12.5  torr  and  mass  flow  < 9.48E-07  lbm/s 

- Thrust  > 2.55  mN  at  Pinlet  > 4.2  torr  and  mass  flow  < 4.85E-06  lbm/s 


• Microthruster  Characterization 

- Thrust 

•^sp 

• Added  Value 


Conical  Nozzle  vs.  Sharp  Edged  Orifice 

Analysis  of  change  in  piston  and  valve  seat  design 
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Mass  Flow  Rate  x 10  A 7 (kg/s) 
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Flow  Direction 


Grid  / Geometry 
for  Flight  Nozzle 
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Mass  Flow  Rate  x 10A4  (grams/sec) 
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Write  detailed  NASA  Tecnical  Memos  or  Technical  Papers 


Subject:  A parametric  study  of  a plug  nozzle,  using  the  Liquid  Propellant  Program  (LPP)  Code 
By:  Stuart  S Dunn,  Douglas  E Coats,  Software  and  Engineering  Associates,  Inc. 


Abstract 

The  Liquid  Propellant  Program  (LPP)  computer  code  is  a super-set  of  the  industry 
standard  Two  Dimensional  Kinetics  (TDK)  computer  code,  which  has  been  developed  by 
Software  and  Engineering  Associates,  Inc.  (SEA,  Inc.)  over  the  past  twelve  years.  The  TDK 
code  uses  a Two-Dimensional  Method  of  Characteristics  solution  with  fully  coupled  finite  rate 
kinetics  for  axially  symmetric  nozzles.  The  chemical  reactions  are  modeled  with  a generalized 
reaction  package  that  includes  3rd  body  efficiencies  and  four  reaction  rate  forms.  The  code 
performs  optional  solutions  for  frozen  or  equilibrium  flow.  TDK  evaluates  discrete  shocks,  both 
attached  or  induced.  The  Transonic  module  models  variable  mixture  ratio  profiles  from  the 
combustion  chamber  injector.  The  Mass  Addition  Boundary  Layer  module  (MABL)  calculates 
the  boundary  parameters  with  the  same  chemistry  options,  and  includes  transpiration  or  tangential 
slot  injection  of  gas  at  the  wall. 

The  LPP  upgrades  include:  planar  nozzles,  scarfed  nozzles,  plug  nozzles,  and  scramjet 
nozzle  configurations.  The  code  evaluates  both  upper  and  lower  wall  flow  simulation,  and 
includes  the  interaction  with  the  external  flow.  The  MABL  module  evaluates  equilibrium 
radiation  heat  transfer  for  both  upper  and  lower  walls.  In  addition,  the  LPP  code  models 
combustion  effects  due  to  injector  inefficiencies  with  the  Spray  Combustion  Analysis  Program 
(SCAP)  module.  The  LPP  package  provides  extensive  post  plotting  capabilities  for  flow 
visualization.  The  LPP  is  sufficiently  fast  and  robust  to  provide  performance  predictions  for 
extensive  parametric  studies  and  sufficiently  accurate  to  provide  flow  field  and  performance 
solutions  for  detailed  studies. 

The  evaluation  of  a planar  or  axially  symmetric  plug  nozzle  has  received  recent  interest 
due  to  the  SSTO  studies.  The  LPP  code  allows  easy  modeling  of  a plug  nozzle  configuration, 
since  the  user  is  allowed  to  input  an  arbitrary  inner  and  outer  wall  geometry  (referred  to  as  the 
plug  and  the  cowl).  The  transonic  analysis  models  both  planar  or  axially  symmetric  annular 
flow,  including  straited  and  variable  mixture  ratio  profiles.  When  the  internal  flow  reaches  the 
exit  of  the  outer  wall,  a Prandtl-Meyer  fan  allows  the  flow  to  expand  to  the  external  pressure. 
At  this  point,  a pressure  boundary  condition  is  applied  for  either  quiescent  sub-sonic,  or  super- 
sonic external  flow.  The  MABL  analyses  is  subsequently  performed  to  evaluate  the  boundary 
layer  losses  for  both  the  inner  and  outer  walls.  Following  JANNAF  standard  procedures,  the 
characteristic  analysis  is  automatically  repeated  with  the  boundary  layer  compensated  wall 
geometry. 

The  above  procedure  was  employed  to  parametrically  evaluate  the  performance  of  several 
plug  nozzle  configurations  at  different  flight  conditions.  The  altitude  compensating  effects  are 
evaluated  and  related  to  ideal  conventional  nozzle  performance.  An  optimization  technique  is 
presented,  which  includes  chemistry,  divergence,  and  boundary  layer  effects.  Graphical  output 
includes  flow  field  contours,  and  wall  property  profiles. 
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SEA  4-26-95/13  th  Propulsion  CFD  Workshop 
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WHAT  IS  LPP? 

The  Liquid  Performance  Program  (LPP)  is 
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WHAT  IS  LPP?  (Continued 
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WHAT  IS  LPP?  (Continued) 
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Flow  Characteristic 
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SSTO  Base  Line  Plug  Noz 
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Isp  Variation  With  Cowl  Length 
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SST0  H2/02  PLUG  NOZZLE  ...  MACH  1 AIR  AT  8 KFT  : PRESSURE 


on  cn  ld 
CD  <nj  — < 
r-  • • 
• o oo 
on  . — . (\J 

<M  CD  O 


CD  CM  CD  CD 
CM  O CD  CO 
CT)  O • • 

■ • —i  CD 

CSJ  OD  CM  LO 

- LD  CD  M 


LD  cn  (XI  LD 
r~  CM  O LD 
CM  CM  • • 

■ • r-  cd 

CM  CO  ' "sf 
O T OO  CM 


cn  on  cn 


OO-Ofr 


00*99 


00*29 


803 


24.00 


804 


MACH  1 AIR  AT  8 KFT 


OCTOt?  00 ' 9£  00  ‘ 2£  00'82  00' VZ 


805 


•00  4.00  8.00  12.00  16.00  20.00  24.00  28- 

X/R* 


STILL  AIR  AT  8 KFT 


00‘  0V  00-90  00-09  00" QZ  00*  7z 


806 


■00  4.00  8-00  12.00  16-00  20-00  24-00  28. 

X/R* 


MACH  1 AIR  AT  8KFT 


00-  Ofr 


00 ’ 9£  00’2£  00-82 


00'fr2 


807 


STILL  AIR  AT  8 KFT 


808 


.00  4.00  8-00  12.00  16-00  20.00  24-00  28. 

X/R* 


Q.  O O 
m O lO 
— 10 


o 

o 

'sf 


o 

LO 

CO 


o 

o 

CO 


o 

LO 

CM 


o 

LO 


809 


o 

o 

CM 


SEA  4-26-95/ 


Isp  Variation  With  Pressure  Rati© 


Q, 

m 


o 

o 

o 

o 

o 

o^“ 

o 

IX) 

o 

IX) 

o 

IT) 

IX) 

CO 

CO 

CM 

810 


CL 

0 

-C 

CO 

Jai 

u. 

1 

Q 

LL. 

o 

c 

o 

cn 

3 

CL 

o 


CO 

25 

a> 

i 

CD 

CM 


55 

c/> 


C/5 

Z 

o 


C/5 

U> 

hJ 

U 

Z 

o 

u 


<2  ^ 


£ 

o 

U 


£ — 

f-H  N 

Oh  O 

b— J L_j 

£ £P  o 

Q r3  ^ 

^ E W 

*•*  ■s 

W • f-H 

0)  ^ 

o ^ 

g 

o3 


X) 

< 


<u 

Z3 

H 


C/5 

H « 

H PQ 


S-H 

o 

T3  4h 

£ S 

g Oh 
-b  a? 

c/3  +3 

C 3 
o pH  _> 

e a ^ 

U o 

Quo; 


C/3 

(D 

O 

§ 

s 

o 


£ 

jD 

E 

0) 

H 60 

aj  Oh 

QJ 


73  ^ 
<U  r2 

> r-1 

o , 

?-H 

P«' 

S g 
M U 

J2  &h 


0) 


<D 

o 

G 

G 


CD 

Oh 


T3 

Jh 

13 

G 

G 

G 

40 

U 

C/3 


3 3 

.2  O 

03  3 

C/3  5-h 

a «2 

<5  ^ 

rv  S3 
cf  D 

S P^ 

O <D 

U Q 

u >, 

T3  *33 

•S.S? 

^ o 


2h  ' ' o 

-b  go  u 

H 'S 

^ o r] 

H £ £ 

T3  bX}.^P 

<U  3 pH 

O M_l  <D 

Z3  ^ -C 

oo  O H 


811 


SEA  4-28-95/13  th  Propulsion  CFD  Workshop 


d) 

CQ 

fi 

o 

4— > 

a, 

O 

Vh 

O 

4-d 

o 

0) 

•c? 

P 

a 

® fH 

o3 


3 

< 

c 

< 

cd 


Oh 

Oh 

hJ 

t2 

73 

d) 

73 

73 

< 


73 

d> 

73 

73 

< 

CD 

PQ 

13 

73 

O 


£ 

o 

E 

d> 
c n 
aj 

PQ 

jD 


00 

< 

03 

^3 

H 


Ph 

Oh 

hJ 

e2 


812 


SEA  4-26-95/13  th  Propulsion  CFD  Workshop 


CFD  ANALYSIS  OF  MODULAR  THRUSTERS  PERFORMANCE 


Ronald  J.  Ungewitter,  James  Beck,  Andrew  Ketchum 
Rocketdyne  Division,  Rockwell  International 
Mail  Code  IB39,  6633  Canoga  Avenue 
Canoga  Park,  CA  91303 

ABSTRACT 

The  effective  performance  of  modular  thrusters  in  an  aerospike  configuration  is 
difficult  to  determine.  Standard  analytical  tools  are  applicable  to  conventional  nozzle 
shapes,  but  are  limited  when  applied  to  an  aerospike  nozzle  (An  aerospike  nozzle  is 
an  altitude  compensating  external  nozzle).  Three  baseline  nozzle  shapes  are  derived 
using  standard  analytical  procedures.  The  baseline  nozzles  sizes  are  restricted  to  fill 
a volume  envelope.  The  three  shapes  are  an  axi-symmetric  round  nozzle,  a 2D  planar 
square  exit  nozzle,  and  a super  elliptic  round  to  nearly  square  nozzle.  The  integrated 
(thruster  /aerospike)  performance  of  the  three  nozzles  is  determined  through  the  use  of 
3-D  viscous  CFD  calculations  where  complex  features  of  the  flowfield  can  be 
accurately  captured.  The  resulting  installed  performance  is  then  used  to  evaluate  the 
efficiency  of  these  nozzle  shapes  for  aerospike  applications. 

The  determination  of  effective  performance  of  a thruster  nozzle  integrated  into 
an  aerospike  nozzle  requires  the  solution  of  the  three  dimensional  turbulent  Navier- 
Stokes  equations.  The  model  used  in  this  study  consisted  of  two  zones;  one  of  the 
upstream  thruster  cowl  surface  so  freestream  conditions  can  be  accurately  predicted, 
and  two,  the  aerospike  surface  beginning  with  with  thruster  outflow  and  extending  to 
the  end  of  the  aerospike  surface.  The  numerical  grid  consisted  of  over  120,000  nodes 
and  used  symmetry  on  the  thruster  centerline  and  edge.  A two  species  non-reacting 
chemistry  model  was  used  to  capture  the  variation  of  fluid  properties  between  the  hot 
plume  gas  and  freestream  air. 

From  the  results  of  the  three  baseline  nozzle  aerospike  calculations,  the 
effective  performance  of  the  nozzle  was  determined.  The  flowfield  of  these 
calculations  do  show  some  variation  between  the  cases.  Recirculation  zones  on  the 
cowl  surface  is  predicted  for  the  2D  planar  nozzle  and  a smaller  one  for  the  super 
elliptic  nozzle.  The  recirculation  is  caused  by  the  strong  pressure  gradient  between 
the  plume  and  freestream  flows.  The  axi-symmetric  nozzle  results  indicates 
recirculation  zones  on  the  thruster  face.  These  recirculation  zones  smooth  the 
pressure  gradient  between  the  plume  and  freestream  flow  limiting  the  formation  of 
recirculation  on  the  cowl  surface.  Thruster  to  thruster  interaction  is  evident  for  the  axi- 
symmetric  and  super  elliptic  calculation  while  the  2D  planar  nozzle  did  not  have  any 
lateral  expansion  in  the  nozzle  so  thruster  to  thruster  interaction  is  limited.  The 
integrated  performance  results,  at  the  altitude  choosen,  shows  very  little  variation 
between  the  three  thruster  shapes.  This  result  allows  for  nozzle  shape  determination 
to  based  on  additional  considerations  (thermal,  structural,  weight)  besides 
performance. 
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BASELINE  NOZZLE  DEFINITIO 
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BASELIME  COMPARISONS 

BASELINE  NOZZLE  AND  ANALYTICAL  METHOD  COMPARISON 
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NOZZLE 

SHAPE 
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2-D  PLANAR 
3-D  SUPER  ELLIPTIC 
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GEOMETRY  OF  INTEGRATED  THRUST  CELL  MODEL 
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THRUST  CELL  TECHNOLOGIES:  INTEGRATED  THRUST  CELL  / AEROSPIKE  ANALYSIS 
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PROPELLANT  CHEMISTRY  FOR  CFD  APPLICATIONS 


,hK, 


R.C.  Farmer,  P.G.  Anderson,  Gary  C.  Cheng  /r  / 


SECA,  Inc. 
Huntsville,  AL 
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ABSTRACT 

Current  concepts  for  Reusable  Launch  Vehicle  design  have  created  renewed  interest  in 
the  use  of  RP-1  fuels  for  high  pressure  and  tri-propellant  propulsion  systems.  Such  designs 
require  the  use  of  analytical  methodology  which  accurately  accounts  for  the  effects  of  real  fluid 
properties,  combustion  of  large  hydrocarbon  fuel  molecules,  and  the  possibility  of  soot 
formation.  These  effects  are  inadequately  treated  in  current  computational  fluid  dynamics  (CFD) 
codes  which  are  used  for  propulsion  system  analyses. 

The  objective  of  this  investigation  is  to  provide  an  accurate  analytical  description  of 
hydrocarbon  combustion  thermodynamics  and  kinetics  which  is  sufficiently  computationally 
efficient  to  be  practical  design  tool  when  used  with  CFD  codes  such  as  the  FDNS  code. 

A rigorous  description  of  real  fluid  properties  for  RP-1  and  its  combustion  products  will 
be  derived  from  the  literature  and  from  experiments  conducted  in  this  investigation.  Upon  the 
establishment  of  such  a description,  the  fluid  description  will  be  simplified  by  using  the 
minimum  of  empiricism  necessary  to  maintain  accurate  combustion  analyses  and  including  such 
empirical  models  into  an  appropriate  CFD  code.  An  additional  benefit  from  this  approach  is  that 
the  real  fluid  properties  analysis  simplifies  the  introduction  of  the  effects  of  droplet  sprays  into 
the  combustion  model. 

Typical  species  compositions  of  RP-1  have  been  identified,  surrogate  fuels  have  been 
established  for  analyses,  and  combustion  and  sooting  reaction  kinetics  models  have  been 
developed.  Methods  for  predicting  the  necessary  real  fluid  properties  have  been  developed  and 
essential  experiments  have  have  been  designed.  Verification  studies  are  in  progress,  and 
preliminary  results  from  these  studies  will  be  presented.  The  approach  has  been  determined  to 
be  feasible,  and  upon  its  completion  the  required  methodology  for  accurate  performance  and  heat 
transfer  CFD  analyses  for  high  pressure,  tripropellant  propulsion  systems  will  be  available. 
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KP-1  COMBUSTION  CHEMISTRY 


RP-1  combustion  model  w/soot  formation 

Verify  model  w/data  from  literature 

Verify  model  w/new  test  data 

Describe  real  fluid  thermodynamic  properties 

Add  real  fluid  HC  properties  to  CICM 

Real  fluid  single  element  model 
Real  fluid  HC,  H2,  and  02 
. Real  fluid  tri-propellant 
Additional  turbulence  models 

Account  for  radiation 


Assemble  elements  to  make  injector  model 


RP-1 


RP-1  is  a straight  run  kerosene  fraction  which  is  subjected  to  acid  washing  and  sulfur 
dioxide  extraction.  The  mean  molecular  weight  is  about  174,  and  the  H/C  ratio  is  about  1.9. 
This  implies  that  RP-1  is  a multicomponent  hydrocarbon  fuel  without  a specified  species 

distribution. 


Composition  by  Hydrocarbon  Typ 

& ,\Y 

Composition 

„ Volume  % ■> . 

dinuclear  aromatics 

2.3 

mononuclear  aromatics 

15.1 

dicyclo-paraffms 

11.6 

monocycyclo-paraffins 

33.8 

branched-paraffms 

15.1 

normal  paraffins 

22.1 

100.0 

or 


Composition 

Volume  % 

aromatics 

17.4 

cyclo-paraffins 

45.4 

paraffins 

37.2 

100.0 
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Distillation  Curve  for  RP-1 
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HYDROCARBONS  ISOLATED  FROM  ONE  REFERENCE  PETROLEUM 


Molecular 

Formula 


Q 


CioHjg 


CioHlg 


C-loH14 


c10h14 


CioHu 


CioHi4 


CnHjo 


Name 

Type 

Normal 
Boiling  Point 

1 .3-Diethylbenzene 

Benzene 

181.10 

l-Methyl-3-propylbenzene 

Benzene 

181.80 

n-Bulylbenzene 

Benzene 

183.27 

Bicycloparaffin 

Bicycloparaffin 

183.4 

Bicycloparaffin 

Bicycloparaffin 

183.7 

1 -Methyl-4-propylbenzene 

Benzene 

183.30 

1,  2-Diethylbenzene 

Benzene 

183.42 

1 ,3-Dimethyl-5-ethylbenzene 

Benzene 

183.58 

1,  4-Diethylbenzene 

Benzene 

183.78 

1,  Methyl-2-propylbenzene 

Benzene 

184.80 

1,  4-Dimethyl-2-ethylbenzene 

Benzene 

186.83 

tr  ans-D  ecahydro-naphthal  ene 

Bicycloparaffin 

187.25 

1 , 3-Dimethyl-4-ethylbenzene 

Benzene 

188.20 

1,  2-Dimethyl-4-ethylbenzene 

Benzene 

| 

Tricyclo  (3. 3. 1.1)-  decane 

Tricycloparaffin 

190 

1,  3-Dimethyl-2-ethylbenzene 

Benzene 

190.01 

1 -Methyl  indan 

Aromatic- 

cycloparaffin 

190.6 

2-Methylindan 

Aromatic- 

cycloparaffin 

191.4 

1,  2-Dimethyl-3  ethylbenzene 

Benzene 

193.91 

cis-Decahydro-naphthalene 

Bicycloparaffin 

195.69 

n-Undecane 

Normal  paraffin 

195.89 

1 ,2,4,5-Tetramethyl-benzene 

Benzene 

196.80 

1 ,2,3 ,5-Tetramethyl -benzene 

Benzene 

198.00 

Bicycloparaffin 

Bicycloparaffin 

202.5 

l-Methyl-3-n-butylbenzene 

Benzene 

204.1 

1,2,3,4-Tetramethyl-benzene 

Benzene 

205.04 

4-Methyl  indan 

Aromatic- 

cycloparaffin 

205.5 

HYDROCARBONS  ISOLATED  FROM  ONE  REFERENCE  PETROLEUM  (Continued) 


C11H16 


C„  H16 


c„h16 


C10Hg 


C„HW 


c„h14 


C 14^30 


C„HU 


Cl3H28 


Mi^io 


C„H10 


C15H32 


C14H30 


C]2H10 


C]3H12 


Cj2H12 


C12Hj2 


Ci2Hi2 


C]2H[2 


Ci2H12 


^ 1^-12 


c12h12 


cJ2h12 


cmhI4 


Name 


1 ,3-Dimethyl -4-n  propylbenzene 


Normal 
Boiling  Point 


Benzene 


1 ,2-Dimethyl-4-n-propylbenzene 


Trimethylethylbenzene 


n-Dodecane 


Naphthalene 


2-Methyl-(l,2,3,4-tetra- 

hydronaphthalene) 


6-Methyl-(l,2,3,4-tetra- 

hydronaphthalene) 


2,6, 10-Trimethylundecane 


5-Methyl-(  1 ,2 , 3 ,4-tetra- 
hydronaphthalene) 


n-Tridecane 


2-Methylnaphthalene 


1-Methylnaphthalene 


2,6,10-Trimethyldodecane 


n-Tetradecane 


Biphenyl 


2-Methylbiphenyl 


2-Ethylnaphthalene 


1 -Ethyl  naphthal  ene 


2 , 6-D  imethyl  naphthal  ene 


2 ,7-D  imethyl  naphthal  ene 


1 ,7-Dimethylnaphthalene 


1 ,6-Dimethylnaphthalene 


1 ,3-Dimethylnaphthalene 


1 ,5-Dimethylnaphthalene 


2,5-Dimethylbiphenyl 


Benzene 


Benzene 


Normal  paraffin 


Dinuclear  aromatic 


Aromatic-cycloparaffin 


Aromatic-cycloparaffin 


Branched  paraffin 


Aromatic-cycloparaffin 


Normal  paraffin 


Dinuclear  aromatic 


Dinuclear  aromatic 


Branched  paraffin 


Normal  paraffin 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


Dinuclear  aromatic 


2 

1,2,3,4-Tetra-hydronaphthalene 

Aromatic-cycloparaffin 

1 

7 

249 


253.52 


25 


25. 


258.7 


18 


ihei 
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RP-1  SURROGATE  FUELS 

Formula 

Species  " • -r; 

Mol  % 

NBPfC) 

Cl3H12 

methylbiphenyl 

17.4 

255 

c12h24 

n-heptylcyclopentane 

45.4 

224 

Ci2H28 

n-tridecane 

37.2 

235 

100.0 

= 173.9 
H/C  = 1.922 


Note:  Critical  pressure  for  RP-1  is  340  psia, 
and  critical  temperature  is  679°K. 
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RP-1  COMBUSTION  PROPERTIES 


HEAT  OF  COMBUSTION 
HOC  = -18640  (Btu/lbm)  = -10.346  (kcal/gm) 


HEAT  OF  VAPORIZATION 
HOY  = 106  (Btu/lbm) 


HEAT  OF  FORMATION  (HOF) 

To  determine  a HOF  for  mixtures,  an  effective  molecular  formula  must  be  specified  and 
used  to  evaluate  the  HOF.  Frequently,  an  arbitrary  molecular  weight  of  100  gms  is  assumed 
for  a basis  for  thermodynamic  calculations. 


H/C 

Source 

Implied  HOC  kcal/g 

1.8624 

-36.01 

Lockheed 

-10.287 

1.90 

-41.6 

SAIC 

-10.305 

2.0 

-44.36 

Aerojet  ’78 

-10.403 

1.9063 

-42.0 

TMX  - 1783 

-10.305 

1.9423 

-33.068(v) 

SP-273 

-10.441(v) 

1.9423 

-38.946(1) 

SP-273 

-10.382(1) 

1.922 

-39.915 

Surrogate  Fuel 

-10.346 
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SPECTROSCOPIC  ANALYSIS 
OF  H-l  GAS  GENERATOR  EXHAUST  GASES 


Cu2 


h2o 


c2h6 


Sample  #1 

Sample  #3 

Sample  #4 

26.15* 

36.82 

37.93 

8.44 

8.86 

9.96 

1.97 

2.22 

2.34 

2.92 

4.20 

4.40 

4.16 

5.62 

4.27 

5.93 

6.93 

7.70 

1.76 


1.58 


0.71 


0.74 


1.14 


1.05 


.99 


5.22 


17.58 


1.73 


0.57 


0.99 


0.97 


1.2 


.3 


2.06 
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Generalized  Combustion  Kinetics  Model 
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Soot  Model 


Initial  Polycyclic  Aromatic  Hydrocarbon  (PAH)  Formation 

• Pyrolysis  and  combustion  of  fuel  to  form  benzene  and  acetylene 

• Implicit  finite-rate  chemistry 

Planar  PAH  Growth 

• HACA 

• Hydrogen  abstraction,  carbon  addition  (through  reactions  with 
acetylene) 

• Oxidation 

• Implicit  finite-rate  chemistry  using  reactions  from  Frenklach,  et  al 

• Properties  for  PAH  compounds 

• Benson’s  group  contribution  method  to  obtain  Cp  for  ideal  gases 

• Benson  data  in  tabular  form  for  300K  < T < 1500K 

• Used  CEC  data  for  selected  species  to  generate  group  contributions 
as  functions  of  temperature  for  300K  < T < 5000K 

• Generated  needed  Cp  data  in  CEC  format 

• S°298  corrected  for  symmetry  and  optical  isomers  using  Benson’s 
data. 


Soot  Model  (Cont.) 
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Cp’s  for  Diphenyl  Compounds 

Dotted  lines  are  synthesized  data 
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Toluene  and  Iso-Octane  Quasiglobal  Kinetics  Model 
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SOOT  PARAMETERS 
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[soot]  = mass  of  soot/volume  of  gas  (g/cm3). 


Elementary  Reactions 


ELEMENTARY  MECHANISM 
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RP-1  COMBUSTION  VERIFICATION 
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MASS  OF  SOOT 


Effect  of  Pressure:  (a)  n-heptane;  (b)  n-dodecane. 
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SOOT  CONCENTRATION,  mg/I  iter  on  dry  basis 


PREDICTION 


A PURE  TOLUENE  DATA  [3-7] 

O 87.5/2.5  BLEND  DATA  [3-7] 


EQUIVALENCE  RATIO 


Comparison  Of  Quasiglobal  Model  Prediction  Of  Soot  Emissions 
With  Experimental  Data  For  Toluene  And  Toluene/Iso-Octane 
Blends.  Computation:  Well-Stirred  Reactor;  Experiment: 
Jet-Stired  Combustor. 
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1.4  1.5  1-6 

EQUIVALENCE  RATIO, 


JSC  Soot  Emissions  For  Toluene/Methylnaphthalene  Blends; 
T = 1900  K;  r = 3 ms. 
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Soot  Concentration  (mg/1  9 300K,  1 atm) 
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Fhotograph  of  the  exit  flow  at  an  exit  pressure  of  2 atm.  Region  1 is  the  undisturbed,  cone, 
region  2 is  the  Prandtl-Msyer  expansion  zone,  and  region  3 is  the  mixing  zone.  The 
spectrometer  line  of  sight  is  indicated  by  a dotted  line. 
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Wavelength  Cum) 


Soot  Pensity  from  GBC  Experiments 
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Rayleigh  Theory  for  Small  Particles 
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Thermo-Kinetics  Characterization  of  Kerosene/RP-1  Combustion 
for  Tripropellant  Engine  Design  Calculations 
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Physical-thermo-chemical  properties  of  kerosene/RP-1. 
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Comparison  of  thermo-chemical  characterization  of  model  fuel  with  reported  data. 
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Mole  Fractions 


Mixture  Ratio  (Oxidizer/Fuel) 


Fig.  1 Calculated  chamber  gas  composition. 
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Mole  Fractions 


Fig.  2 Calculated  nozzle  exit  gas  composition. 


872 


Temperature,  deg  K 
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Kerosene/RP- 1 quasi-global  combustion  kinetics  mechanism.  Kf  = AT  e 
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M stands  for  third-body  collision  partner 


Kerosene  Substitute  Fuel  Model/Quasi-gioba! 
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Multiphase  FDNS  Numerical  Methodology  and  Physical  Models 

for  Triaxial  Tripropellant  Injectors 
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- real  gas:  kerosene  and  8 other  species 

- real  liquid:  LOX  (NBS  table) 

: kerosene  (Lefebvre) 
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Application  of  Optimization  Techniques  to  Design 
©f  Unconventional  Rocket  Nozzle  Configurations  / a - )//%/ 

W.  Follett,  A.  Ketch  urn,  A.  Darian,  Y.  Hsu  /Of 

Rocketdyne  Division,  Rockwell  International,  Canoga  Park,  California 

Abstract 

Several  current  rocket  engine  concepts  such  as  the  bell-annular  tripropellant  engine,  and  the 
linear  aerospike  being  proposed  for  the  X-33,  require  unconventional  three-dimensional 
rocket  nozzles  which  must  conform  to  rectangular  or  sector  shaped  envelopes  to  meet 
integration  constraints.  These  types  of  nozzles  exist  outside  the  current  experience  database, 
therefore,  application  of  efficient  design  methods  for  these  propulsion  concepts  is  critical  to 
the  success  of  launch  vehicle  programs. 

The  objective  of  this  work  is  to  optimize  several  different  nozzle  configurations,  including  2-0 
and  3-D  geometries.  Methodology  includes  coupling  CFD  analysis  to  genetic  algorithms  and 
Taguchi  methods,  as  well  as  implementation  of  a streamline  tracing  technique.  Results  of 
applications  are  shown  for  several  geometry  classes  including:  3-D  thruster  nozzles  with 
round  or  superelliptic  throats  and  rectangular  exits,  2-D  and  3-D  thrusters  installed  within  a 
bell  nozzle,  and  3-D  thrusters  with  round  throats  and  sector  shaped  exits. 

Due  to  the  novel  designs  considered  for  this  study,  there  is  little  experience  base  which  can 
be  used  to  guide  the  effort  and  limit  the  design  space.  With  a nearly  infinite  parameter  space 
to  explore,  simple  parametric  design  studies  cannot  possibly  search  the  entire  design  space 
within  the  time  frame  required  to  impact  the  design  cycle.  For  this  reason,  robust  and  efficient 
optimization  methods  are  required  to  explore  and  exploit  the  design  space  to  achieve  high 
performance  engine  designs.  Five  case  studies  which  examine  the  applications  of  various 
techniques  in  the  engineering  environment  are  presented  in  this  paper. 


The  first  study  uses  two-dimensional  CFD  coupled  to  Taguchi  methods  to  determine  optima! 
design  parameters  for  the  D-1  test  engine  being  built  for  the  SSTO  Advanced  Propulsion 
Technology  contract.  The  D-1  engine  utilizes  a ring  of  small  thrusters  within  a larger  bell 
nozzle.  This  study  was  used  to  determine  the  optimal  value  of  four  design  variables  to 
achieve  the  best  overall  performance  during  both  low  altitude  (thrusters  firing)  and  high 
altitude  (thrusters  not  firing)  operational  modes.  Two  other  case  studies  investigate  the 
problem  of  using  multidisciplinary  techniques  to  optimize  a 3-D  thruster  design  with  both 
genetic  algorithms  and  Taguchi  methods.  The  relative  strengths  and  weaknesses  of  these 
two  methods  are  apparent  when  using  them  to  solve  this  problem  using  up  to  21  design 
variables.  This  thruster  is  also  designed  using  streamline  tracing  techniques  for  the  fourth 
case  study. 

The  final  study  uses  Taguchi  methods  to  determine  the  optimal  3-D  thruster  module  design 
when  installed  in  a bell  nozzle.  This  requires  full  3-D  solutions  of  the  thruster  and  bell  nozzles 
to  quantify  module-to-module  interaction  effects. 

Software  which  couples  optimization  techniques  to  CFD  have  tremendous  potential  as 
aerodynamic  design  tools.  However,  to  function  effectively  in  the  engineering  environment, 
the  optimization  algorithms  must  be  robust  and  efficent.  Several  optimization  techniques 
have  been  demonstrated  for  rocket  nozzle  design,  and  their  performance  on  these  real  world 
applications  has  been  assessed. 
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3-D  THRUST  CELL  OPTIMiZATIO 
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3-D  STREAMLINE  TRACING  METHOD 
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3-D  BELL-ANNULAR  D-1  INSTALLED  THRUSTER  DESIGN 
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Transonic  Aerodynamic  Characteristics  of  A Proposed 
Wing  Body  Reusable  Launch  Vehicle  Concept 

Anthony  M.  Springer 

NASA  Marshall  Space  Flight  Center  ■ ' c 

/ L3  cD  / o ' 

A proposed  wing  body  reusable  launch  vehicle  was  tested  in  the  NASA  Marshall  Space  Flight  / /*k'5’  j 

Center  14X14-inch  trisonic  wind  tunnel  during  the  winter  of  1994.  This  test  resulted  in  the 
vehicle's  subsonic  and  transonic,  Mach  0.3  to  1 .96,  longitudinal  and  lateral  aerodynamic 
characteristics.  The  effects  of  control  surface  deflections  on  the  basic  vehicles  aerodynamics 
including  a body  flap,  elevons,  ailerons,  and  tip  fins  are  presented. 

As  an  outcome  of  NASA’s  1993  Access  to  Space  study,  a more  in-depth  follow  on  study  was 
undertaken.  Three  candidate  reusable  launch  vehicle  configurations  which  would  provide 
reusable  single  stage  to  orbit  capability  were  selected.  A wing  body  configuration  was  one  of 
these  candidate  concepts,  the  other  two  concepts  being  a vertical  lander  and  a lifting  body.  The 
wing  body  configuration  was  a direct  outgrowth  of  the  access  to  space  option  three  reference 
single  stage  to  orbit  rocket  vehicle.  This  vehicle  matured  during  the  subsequent  reusable  launch 
vehicle  (RLV)  study  into  the  vehicle  which  was  tested.  Initially,  the  vehicles  aerodynamic 
characteristics  were  determined  using  aerodynamic  prediction  codes.  To  obtain  a better  fidelity 
in  the  aerodynamic  data,  a series  of  scale  models  of  the  proposed  wing  body  vehicle  were  tested 
at  the  NASA  Marshall  Space  Flight  Center  (MSFC)  and  the  NASA  Langley  Research  Center 
(LaRC).  The  vehicle  was  tested  at  low  subsonic  and  hypersonic  conditions  at  LaRC  and  at 
subsonic,  transonic,  and  supersonic  conditions  at  MSFC.  The  results  of  the  transonic  testing  in 
MSFC's  14-Inch  Trisonic  Wind  Tunnel  (TWT)  facility  are  presented  herein. 

A .004  scale  RLV  wind  tunnel  model  was  tested  during  the  winter  of  1994  at  the  NASA  Marshall 
Space  Flight  Center  14X1 4-inch  trisonic  wind  tunnel  (TWT).  The  subsonic  and  transonic, 

Mach  0.3  to  Mach  2.0,  aerodynamic  characteristics  of  the  WB001  reference  wing-body  vehicle 
were  determined.  This  wind  tunnel  test  provided  aerodynamic  data  for  the  basic  vehicle,  wing 
and  body  contributions,  and  control  surface  increments.  The  data  derived  from  this  test  were 
used  to  construct  an  aerodynamic  database  for  flight  mechanics  and  structural  loads  studies  on 
the  wing  body  vehicle. 

The  WB001  vehicle  is  generically  a wing-body  combination.  The  body  consists  of  a drooped  nose 
followed  by  a cylindrical  core  section  28.55  ft  in  diameter,  full  scale,  with  a total  body  length 
of  185.6  ft,  full  scaie.  The  wing  is  a NACA-0010  airfoil  at  the  root  linearly  varying  to  a NACA- 
0012  airfoil  at  the  tip  with  a 54  degree  leading  edge  sweep,  3.5  degrees  of  dihedrai,  and  an 
aspect  ratio  of  1.91.  Control  surfaces  for  this  configuration  consist  of  ailerons,  elevons  and  tip 
fins. 

The  vehicle  is  longitudinally  stable  and  can  be  trimmed  at  both  subsonic  and  transonic  Mach 
numbers.  This  assumes  a vehicle  center  of  gravity  at  68.6%  body  length  or  127.32  feet  aft  of 
the  nose.  At  subsonic  Mach  numbers,  the  vehicle  is  stable  in  trim  for  all  control  deflections. 

The  vehicle  for  the  subsonic  Mach  range  can  be  trimmed  at  the  desired  angle-of-attack  for 
entry,  approximately  15  degrees.  This  trim  angle  is  accomplished  through  various  control 
surface  deflections,  see  figure  84.  The  vehicle  for  the  transonic  Mach  range,  Mach  0.95  to  2.0, 
has  stable  trim  points  but  not  at  the  desired  angle-of-attack,  approximately  15  degrees  angle- 
of-attack.  It  can  be  extrapolated  from  the  current  data  trends  that  for  a larger  elevon  deflection 
between  20  and  30  degrees,  the  vehicle  will  be  neutrally  stable  at  the  desired  trim  point  of  15 
degrees. 

The  WB001  vehicle  is  laterally  unstable  for  the  subsonic  and  transonic  Mach  range.  The  tip  fin 
deflections  provide  a trim  angle  range  of  approximately  1 to  2 degrees,  therefore,  larger  tip 
fins  and  deflectable  surfaces  are  desirable.  Enlarging  the  tip  fins  by  an  approximate  geometric 
factor  of  3 to  4 should  result  in  the  vehicle  being  neutrally  stable. 
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Ascent  Aerodynamic  Pressure  Distributions  on  WB001 

B.  Vu,  J.  Ruf,  F.  Canabal  <3^  -P 

CFD  Branch 


J.  Brunty 

System  Load  Branch 

To  support  the  reusable  launch  vehicle  concept  study,  the  aerodynamic  data 
and  surface  pressure  for  WB001  were  predicted  using  three  CFD  codes  at 
several  flow  conditions  during  the  ascent  phase.  The  results  have  been 
compared  between  code  to  code  and  code  to  aerodynamic  database  as  well  as 
available  experimental  data.  A set  of  particular  solutions  have  been  selected 
and  recommended  for  use  in  preliminary  conceptual  designs.  These  CFD 
results  have  also  been  provided  to  the  structure  group  for  wing  loading 
analyses. 
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finite-difference,  pressure-based 

used  at  MSFC  for  reacting  flow  analyses 


- predicted  stagnation  Cp  agrees  with  isentropic  theory 

- M = 1 .1  Cp-1 .36  (CFD  / GASP) 

- M = 1.1  Cp=1.34  (Theory) 
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Assessment  of  Lifting  Body  Linear  Aerospike  Plume  Effects 
on  Vehicle  Aerodynamics 


Joseph  H Ruf  MSFC/ED32 
Alonzo  L.  Frost  MSFC/ED34 
Bruce  Vu  MSFC/ED32 
Francisco  Canabal  MSFC/ED32 


The  lifting  body/linear  aerospike  is  one  of  three  configurations  being  studied 
for  an  SSTO  vehicle.  A preliminary  aerodynamic  database  existed  for  then  current 
lifting  body  configuration,  however,  this  data  base  was  developed  without 
considering  plume  effects.  A combined  effort  by  the  Computational  Fluid  Dynamics 
and  Experimental  Fluids  Dynamics  Branches  was  undertaken  to  determine  first 
order  effects  of  plume /external  flow  interactions  on  vehicle  aerodynamics  of  this 
lifting  body/linear  aerospike  configuration.  Of  interest  were  plume 
pumping/ entrainment  at  low  Mach  numbers  and  plume  induced  separation  of  flow 
over  the  vehicle  at  higher  altitudes.  The  CFD  analysis  included  combinations  of 
four  Mach  numbers,  two  angles  of  attack  and  four  throttle  settings.  The  majority  of 
the  CFD  was  two  dimensional  centerline  analysis  of  the  lifting  body  /aerospike. 
Incremental  plume  effects  were  derived  by  comparing  the  power-on,  power-off,  and 
throttled  cases  and  were  extrapolated  to  the  preliminary  aerodynamic  database. 

The  plume  had  little  effect  on  the  vehicle  aerodynamics  for  supersonic 
freestream  velocities.  At  subsonic  freestream  velocities,  the  plume  affected  the 
vehicle  aerodynamics  through  both  jet  pumping/ entrainment  and  the  jet  flap  effect. 
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CONFIGURATION  K-10 


Sref  = 5,600  sq.  ft. 

Swet  = 13,100  sq.  ft. 

Total  Internal  Volume  = 67,000  cu.  ft. 
Main  Hyd.  Tank  Volume  = 47,500  cu.  ft. 
Fwd  LOX  Tank  Volume  = 5,200  cu.  ft. 
Aft  LOX  Tank  Volume  = 1 2,700  cu.  ft. 
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6 degree  a power-on  vs.  power-off 

» significant  plume  effects  through  plume  entrainment  and  jet  flap  effect 


nach  NuMber 


Mach  5.D  FreestreaM,  126Kft 
KID  2D  Centerline,  0 degrees^-fi 


CONTOUR  LEUELS 


Mach  Number 


K10  2D  Centerline  Pressure  Coefficients 
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K10  2D  Centerline  Pressure  Coefficient 


Mach  0.6  Alpha=0  degrees 


power-off  vs.  power-on 
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K10  2D  Centerline  Pressure  Coefficients 
Mach  0.6  Alpha=0  degrees  Mach  0.6  Alpha=6  degrees 

P-On  vs.  P-On-Throttled  (70%/1 30%  RPL  Upper/Lower  Thrusters)  P-On  vs.  P-On-Throttled  (70%/1 30%  RPL  Upper/Lower  Thrusters) 
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DTC  - differential  throttle  control;  upper/lower  @ 1 30%/70%;  Pon  = 1 00%/1 00% 
dCp  = difference  in  pressure  coefficient  {power  on  - power  off) 


Sref=5600  sq  ft;  Lref=1 450.0  in.;  mom.  ref.  @ 72%. 


Plume  effects  on  LB  pressure  distribution  @ Mach  0.6,  angle-of-attack=6  deg 
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DIG  - differential  throttle  control;  upper/lower  @ 130%/70%;  Ron  = 1 00%/1 00% 
dCp  = difference  in  pressure  coefficient  (power  on  - power  off) 


Plume  effects  on  LB  normal  force  coeff  (M=0.6) 


angle-of-attack  (deg) 

Sref=5600  sq  ft;  Lref=1 450.0  in.;  mom.  ref.  <§>  72%. 
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A methodology  has  been  developed  to  generate  first  order  plume 
effect  increments  for  a power-off  aerodynamic  database  using  2D 
centerline  CFD  analysis. 


TPS  Sizing  for 
Access-to-Space  Vehicles 


William  Henline,  David  Olynick  and  Grant  Palmer 
NASA  Ames  Research  Center,  MS  230-2, 

Moffett  Field,  CA  94035-1000 
Y.-K.  Chen 

Eloret  Institute,  MS  234-1,  Moffett  Field,  CA  94035-1000 


Abstract 

A study  was  carried  out  to  identify,  develop,  and  benchmark  simulation 
techniques  needed  for  optimum  TPS  material  selection  and  sizing  for  reusable 
launch  vehicles.  Fully  viscous,  chemically  reacting,  Navier-Stokes  flow  solu- 
tions over  the  Langley  wing-body  single  stage  to  orbit  (SSTO)  configuration 
were  generated  and  coupled  with  an  in-depth  conduction  code.  Results  from 
the  study  provide  detailed  thermal  protection  system  (TPS)  heat  shield  ma- 
terials selection  and  thickness  sizing  for  the  wing-body  SSTO.  These  results 
are  the  first  ever  achieved  through  the  use  of  a complete,  trajectory  based 
hypersonic,  Navier-Stokes  solution  database.  TPS  designs  were  obtained 
for  both  laminar  and  turbulent  entry  trajectories  using  the  Access-to-Space 
baseline  materials  such  as  tailorable  advanced  blanket  insulation  (TABI). 
The  TPS  design  effects  (material  selection  and  thicknesses)  of  coupling  ma- 
terial characteristics  to  the  aerothermal  enviroment  are  illustrated.  Finally, 
a sample  validation  case  using  the  shuttle  flight  data  base  is  included. 

For  the  laminar  trajectory,  the  TPS  areal  mass  density  is  1.2  lbm/ft2, 
while  the  turbulent  trajectory  yields  slightly  less  than  1.3  lbm/ft2.  An  addi- 
tional conclusion  from  this  study  is  that  the  TABI  blankets  will  have  to  be 
manufactured  in  thicknesses  greater  than  1. 5-2.0  inches.  Further,  if  typical 
turbulent  flow  conditions  are  found  on  these  SSTO  vehicles  during  re-entry, 
some  of  the  baseline  materials  may  experience  significant  over-temperatures. 
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NASA  Ames  Research  Center 
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(Including  operations)  of  all  candidate  TPS  and  thus  OMB  TPS  criteria 


Fully  Coupled  Thermal  Analysis  for  TPS  Sizing 
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Trajectory  (Altitude  - Velocity)  Plot  for  the  LaRC  SSTO  Vehicle 
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fl  ajectoi y Time  (sec)  LaRC  - Weilmcunstcr/Wurslcr  4/94 


Stagnation  Point  Heating  Plot  for  LaRC  SSTO  Trajectory 


969 


(0 

BBBOa 

m 

-8-fl 

CO 

2 

U) 


m 

□ 

hm 

3 

O 

c 

o 

o 

22 

5 

2 

a> 

a 

E 

G) 

S— 

Q) 

O 

re 

t 

</) 


970 


Winged  Body  Configuration  1300  s 
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Winged  Body  Configuration  1300  S 
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Surface  TPS  Thickness  (in.) 
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Top  Layer  TPS  Thickness  (in.)  for  the  LaRC  Winged  Body  SSTO  Vehicle  (Total  Heating  Time,  6200  sec) 


(TURBULENT  FLOW  SOLUTION) 


LaRC  SSTO  Vehicle  Centerline  TPS  Thicknesses 
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Effect  of  TPS  Material  Properties  on  Surface  Temperatures 
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COMPUTATIONAL  ISSUES  ASSOCIATED  WITH  TEMPORALLY 
DEFORMING  GEOMETRIES  SUCH  AS  THRUST  VECTORING 

NOZZLES 


Kishore  Boyalakuntla,  Bharat  K.  Soni,  Hugh  J.  Thornburg  and 

Robert  Yu 

NSF  Engineering  Research  Center  for 
Computational  Field  Simulation 
Mississippi  State  University 
Mississippi  State,  MS  39762 

ABSTRACT 

During  the  past  decade  Computational  Simulation  of  fluid  flow  around  complex  configurations  has  progressed  significantly  and  many 
notable  successes  have  been  reported,  however  unsteady  time-dependent  solutions  are  not  easily  obtainable.  The  present  effort  involues 
unsteady  time  dependent  simulation  of  temporally  deforming  geometries.  Grid  generation  for  a complex  configuration  can  be  a time  con- 
suming process  and  temporally  varying  geometries  necessitate  the  regeneration  of  such  a grid  for  every  time  step.  Traditional  grid  genera- 
tion techniques  have  been  tried  and  demonstrated  to  be  inadequate  to  such  simulations.  NURBS  based  techniques  provide  a compact  and 
accurate  representation  of  the  geometry.  This  definition  can  be  coupled  with  a distribution  mesh  for  a user  defined  spacing.  The  present 
method  greatly  reduces  cpu  requirements  for  time  dependent  remeshing,  facilitating  the  simulation  of  more  complex  unsteady  problems. 
A Thrust  Vectoring  Nozzle  has  been  chosen  to  demonstrate  the  capability  as  it  is  of  current  interest  in  aerospace  industry  for  better  man- 
oeverability  of  fighter  aircraft  in  close  combat  and  in  post  stall  regimes.  This  current  effort  is  the  first  step  towards  multidisciplinary  design 
optimization  which  involues  coupling  the  aerodynamic  heat  transfer  and  structural  anlysis  techniques.  Applications  include  simulation 
of  temporally  deforming  bodies  and  areoelastic  problems. 

A NURBS  based  volume  grid  generation  technique  is  used  for  remeshing  at  each  timestep.  Remeshing  is  easily  accomplished  by  varying 
the  control  points  and  time  dependent  motion  is  contained  in  the  motion  of  the  control  points.  Timestep  controls  the  movement  of  control 
points.  Great  flexibility  in  geometric  definition  is  achieved.  The  grid  generation  code  is  succesfully  coupled  with  UBIFLOW  and  INS3d 
which  are  compressible  and  incompressible  flow  solvers  respectively. 

Various  geometries  such  as  converging  diverging  nozzle,  duct  and  thrust  vectoring  nozzle  have  been  simulated  and  wiH  be  presented. 
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7,  Flowsolver  that  can  handle  moving  boundaries . 
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PRETTY  PICTURES  IN  BETWEEN 
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for  temporally  varying  geometries. 


CONTROL  POINTS  FOR  AIRCRAFT  WING  (4x2x9) 
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Distribution  Mesh 
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wing  generated  with  control  points  and  distribution  mesh 
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operating  on  aircraft  carriers. 


Grid  at  time  steps  of  0,  273,  363,  423(1/2  cycle) 


Pressure  Distribution  for  timesteps  of  0,  273,363,423 
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Circular  to  Rectangular  Transition  Duct 
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OJLAR  TO  RECTANGULAR  TRANSITION  DUCT  ( 0 timestep)  rhrust  VecTonnS  ( 70  timestep) 
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THRUST  VECTORING  ( 150  timesteps) 


0.009  0.096  0.192  0.288  0.383  0.479 


1004 


— -Jz 

Cis  E J 
IUO  52 

y g til 
?« 
o S5  sS 
ZSOH 

Ulc  OE 


C/2 

C-> 

a>  5* 

C/2  QJ 

tH  5/2 
. c/2 

2 8 
® 2 
P A 
Oh  N 

U X 

hH  g 


o 

fl 

s 

o> 

d 

53 

d 

cc 

« 

OD 

V 

d 

s 

n 

d 

GC 

© 

d 

fH 

o 

C/2 

• PM 

-« 

o 

'O 

d 

o 

cc 

o 

a> 

T3 

•PM 

o 

T3 

5* 

OD 

OD 

S3 

©PM 

► 

> 

O 

o 

C/2 

a 

£ 

0> 

o 

1005 


nozzle  have  been  snccesfully  simulated. 


1007 


HYBRID  GRID  TECHNIQUES  FOR  PROPULSION  APPLICATIONS 


Roy  P.  Koomullil,  Bharat  K.  Soni,  and  Hugh  J.  Thornburg 
NSF  Engineering  Research  Center  for 
Computational  Field  Simulation 
Mississippi  State  University 
Mississippi  State,  MS  39762 

ABSTRACT 

During  the  past  decade  computational  simulation  of  fluid  flow  for  propulsion  applications  has  progressed  significantly,  and  many  notable  successes 
have  been  reported  in  the  literature.  However,  the  generation  of  a high  quality  mesh  for  such  problems  has  often  been  reported  as  a pacing  item.  Hence, 
much  effort  has  been  expended  to  speed  this  portion  of  the  simulation  process.  Several  approaches  have  evolved  for  grid  generation.  Two  of  the  most 
common  are  structured  multi-block,  and  unstructured  based  procedures.  Structured  grids  tend  to  be  computationally  efficient,  and  high  aspect  ratio 
cells  necessary  for  efficiently  resolving  viscous  layers.  Structured  multi-block  grids  may  or  may  not  exhibit  grid  line  continuity  across  the  block  inter- 
face. This  relaxation  of  the  continuity  constraint  at  the  interface  is  intended  to  ease  the  grid  generation  process,  which  is  still  time  consuming.  Flow 
solvers  supporting  non-contiguous  interfaces  require  specialized  interpolation  procedures  which  may  not  ensure  conservation  at  the  interface.  Un- 
structured or  generalized  indexing  data  structures  offer  greater  flexibility,  but  require  explicit  connectivity  information  and  are  not  easy  to  generate 
for  three-dimensional  configurations.  In  addition  unstructured  mesh  based  schemes  tend  to  be  less  efficient  and  it  is  difficult  to  resolve  viscous  layers. 
Recently,  hybrid  or  generalized  element  solution  and  grid  generation  techniques  have  been  developed  with  the  objective  of  combing  the  attractive 
features  of  both  structured  and  unstructured  techniques.  In  the  present  work  recently  developed  procedures  for  hybrid  grid  generation  and  flow  simula- 
tion are  critically  evaluated,  and  compared  to  existing  structured  and  unstructured  procedures  in  terms  of  accuracy  and  computational  requirements. 

In  the  present  grid  generation  procedure  multi-body  configurations  are  decomposed  into  a number  of  simple  geometric  entities.  A structured  grid  gener- 
ator is  first  employed  to  construct  a high  quality  grid  around  the  body  with  appropriate  packing.  One  grid  must  be  designated  as  a main  grid  and  enclose 
the  solid  surfaces  of  all  other  component  grids.  Upon  completion  these  structured  grids  are  converted  to  the  hybrid  grid  data  structure  format.  Based 
upon  an  input  normal  distance  from  the  surface,  holes  are  cut  in  the  main  grid  for  each  component  grid.  Overlapping  and  hole  cells  are  deleted  from 
the  hybrid  grid  data  structure.  Delaunay  triangulation  is  then  used  to  construct  cells  to  fill  the  void  between  the  cut  main  grid  and  the  truncated  compo- 
nent grid.  Upon  completion  of  this  procedure  the  hybrid  grid  is  written  in  a format  useable  by  the  flow  solver. 

The  non-dimensionalized  Euler  equations  in  integral  form  provide  the  mathematical  formulation  for  this  scheme.  The  discretized  flow  domain  is  repre- 
sented by  a set  of  non  overlapping  polygons  and  the  cell  averaged  variables  are  stored  at  each  cell  center.  Each  individual  cell  is  treated  as  its  own 
control  volume.  The  numerical  flux  at  the  cell  edge  is  calculated  using  Roe’s  approximate  Riemann  solver.  An  assumed  linear  distribution  in  each 
cell  is  employed  to  reconstruct  the  edge  values,  which  results  in  a second  order  discretization.  The  flux  limiting  procedure  of  Barth  is  used  to  suppress 
spurious  oscillations  near  discontinuities.  An  implicit  pseudo-time  integration  procedure  using  the  Generalized  Minimum  RESidual  (GMRES)  meth- 
od for  solving  the  sparse  matrix  system  is  employed.  The  results  has  been  varified  with  the  standard  benchmark  results. 
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Scramjet  Inlet  Like  Geometry  (NPARC) 
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0.654  2.734  4.815  6.895 
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Scramjet  Inlet  Like  Geometry  ( HYBRID  ) 
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Iteration  Number 
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A Structured  Grid  Based  Solution-adaptive  Technique  for  Complex  Separated 

Flows  / 

by 


Hugh  Thornburg,  Bharat  K.  Soni,  Boyalakuntla  Kishore,  and  Robert  Yu 
NSF  Engineering  Research  Center  for 
Computational  Field  Simulation 
Mississippi  State  University 
Mississippi  State,  MS  39762 


ABSTRACT 


The  objective  of  this  work  has  been  to  enhance  the  predictive  capability  of  widely  used  CFD  codes 
through  the  use  of  solution  adaptive  gridding.  Most  problems  of  engineering  interest  involve  multi- 
block grids  and  widely  disparate  length  scales.  Hence,  it  is  desirable  that  the  adaptive  grid  feature 
detection  algorithm  be  developed  to  recognize  flow  structures  of  different  type  as  well  as  differing 
intensity,  and  adequately  address  scaling  and  normalization  across  blocks.  In  order  to  study  the  ac- 
curacy and  efficiency  improvements  due  to  the  grid  adaptation,  it  is  necessary  to  quantify  grid  size 
and  distribution  requirements  as  well  as  computational  times  of  non-adapted  solutions.  Flowfields 
about  launch  vehicles  of  practical  interest  often  involve  supersonic  freestream  conditions  at  angle 
of  attack  exhibiting  large  scale  separated  vortical  flow,  vortex-vortex  and  vortex-surface  interac- 
tions, separated  shear  layers  and  multiple  shocks  of  different  intensity.  In  this  work  a weight  func- 
tion and  an  associated  mesh  redistribution  procedure  is  presented  which  detects  and  resolves  these 
features  without  user  intervention.  Particular  emphasis  has  been  placed  upon  accurate  resolution 
of  expansion  regions  and  boundary  layers. 

Flow  past  a wedge  at  Mach  = 2.0  is  used  to  illustrate  the  enhanced  detection  capabilities  of  this  newly 
developed  weight  function.  Figure  1 presents  weight  functions  evaluated  using  the  previous  proce- 
dure, lower  half  plane,  as  well  as  the  current  procedure,  upper  half  plane. 


WEIGHT  WITH/WITHOUT  BOOLEAN 


Figure  1.  Comparison  of  Weight  Functions. 


ADAPTED  GRID  WITH/WlTHOUT  (BOOLEAN) 


(80,50,5) 


Figure  2.  Comparison  of  Adapted  Grids. 


It  can  be  observed  that  both  weight  functions  clearly  detected  the  primary  shock.  It  can  also  be  seen 
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that  the  expansion  fan,  boundary  layer,  and  the  reflected  shocks  are  much  more  clearly  represented 
in  the  current  weight  function.  Adapted  grids  using  both  weight  function  formulations  are  presented 
in  Fig.  2.  The  high  gradient  regions  of  the  expansion  region  are  only  reflected  in  the  adapted  grid 
using  the  new  weight  function.  The  reflected  shock  is  also  much  sharper.  Figure  3 compares  the 
solution  obtained  using  the  current  adaption  procedure  with  that  obtained  using  the  original  grid. 
The  enhanced  resolution  is  clearly  evident. 


ADAPTED/ORIGINAL  SOLUTION  (BOOLEAN) 


Supersonic  flow  at  Mach=1.45  and  14  degree  angle  of  attack  has  been  simulated  around  a tangent- 
ogive  cylinder.  The  grid  and  associated  flow  solution  constructed  after  two  adaption  cycles  using 
hybrid  differencing  of  the  grid  equations  and  the  current  weight  functions  is  presented  is  presented 

in  Figure  4. 


Figure  4.  Adapted  grid  after  two  cycles.  Figure  5.  Adapted  grid  after  two  cycles. 


Figure  5 presents  the  grid  constructed  using  the  previous  weight  function  and  the  same  flow  condi- 
tions and  number  of  adaptation  cycles.  Figures  5 and  6 present  streamwise  cuts  of  the  two  grids 
shown  in  Figs  4 and  5 at  X/D  = 5.5  and  7.5  respectively 
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Figure  6.  X/D  = 5.5 


Figure  7.  X/D  = 8.5. 


Figure  8 present  the  flow  solution  obtained  using  the  NPARC  [NASA  1993]  flow  solver,  the  KE 
turbulence  model  option  and  two  adaptation  cycles.  Figure  9 presents  the  associated  weight  func- 
tion. 


Figure  8 Normalized  Stagnation  Pressure.  Figure  9.  Weight  Function. 


Examples  will  presented  to  demonstrate  the  capability  for  solution-adaptive  regridding  of  multi- 
block launch  vehicle  simulations. 
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COMPLEX  GEOMETRY /COMPLEX  PHYSICS 
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OBJECTIVES 


Improved  resolution  of  complex  flows  through  the  use  of  solti 
tion  adaptive  gridding 

1.  Develop  a weight  function  suitable  for  use  with  a solution  adaptive 
grid  redistribution  procedure  for  complex  flows,  including  viscous 
dominated  separation. 

2.  Minimum  user  inputs. 

3.  Appropriate  feature  detection  for  a wide  range  of  flow  features  (Vorti- 
cies,  Shear  layers,  Shocks). 

4.  Robust  redistribution  procedure  for  use  with  weight  function. 


GOVERNING  EQUATIONS  FOR  GRID  MOVEMENT 


1.  Inverted  form: 

Z Z 8ij  + Z Skk  pk  rv  = o 

i—\ j—\  jfc= 1 

Where: 

rjj  : Position  vector, 

g)j  : Contravariant  metric  tensor, 

I'  : Curvilinear  coordinate,  and 
Pk  : Control  Function. 

2.  Control  of  distribution  and  characteristics  of  grid  system  can  be 
achieved  by  varying  control  Functions  Pk- 


W\ 
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6.  For  Multiple  dimensions: 

Pk=Wtk/W,  k = 1,2,3 


CHARACTERISTICS  OF  WEIGHT 

FUNCTIONS 

r 

1 . Weight  functions  approximately  equidistributed 
over  solution  domain. 

2.  Determine  grid  spacing  and  characteristics. 

3.  Approximation  to  local  truncation  error. 

• Use  lower  order  derivatives  to  approximate  high  order  trunca- 
tion error  terms. 

• Detect  structures  of  disparate  strength. 

• Minimum  variation  of  coefficients. 


EVALUATION  OF  WEIGHT  FUNCTIONS 
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maximum  in  order  to  maintain  the  relative  strength. 


ADAPTIVE  GRID  PROCEDURE 


1. Read  PLOT3D  grid  and  solution  files. 

2.  Evaluate  weight  function, 

(no  input  parameters). 

3.  Evaluate  and  smooth  P^. 

4. Integrate  grid. 

5. Interpolate  P^  onto  current  adapted  grid. 

6.  Repeat  steps  4 and  5 until  convergence. 

7.  Ouput  adapted  grid. 


SOLUTION  OF  GRID  EQUATIONS 


1.  Solution  difficulties  transferred  form  flow  equations  to 
grid  equations. 

2.  Accuracy  not  as  important  for  postulated  law. 

3.  Adaptive  Central/Upwind  differencing  scheme,  based 
upon  forcing  function  gradients. 

4.  Integrated  in  time  using  CSIP. 

5.  Non-linear  terms  are  quasi-linearized. 

6.  Explicit  boundary  point  movement. 

7.  Precise  geometry  definition  is  critical. 


Grid  Central  I )i 


T 


Very  important. 

• Orthogonality. 

• Skewness. 

Non  Uniform  Rational  B-Spline  (NURBS) 
representaion  of  arbitrary  surface(Yu). 

Boundary  surface  redistribution  based  on 
specified  region  of  surface. 

• Explicit. 

• Local  iteration  for  desired  distribution. 

• Can  be  used  to  keep  sharp  comers,  and  to  transfer  information 
between  blocks. 
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Figure  2.  Comparison  of  Adapted  Grids. 


Case  1 NPARC  Si 
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Case  1 NP ARC  Grid  7 a2  KE  Model 


Wm 


U BIFLOW  Density,  -light  Marked  Cells 
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SUMMARY 


1.  Developed  Weight  function  which  requires  no 
user  input. 

2.  Implemented  adaptive  upwind/central 
difference  scheme. 

3.  Demonstrated  enhanced  grid  resolution. 

• Thinner  shocks. 

• Stronger  circular  vorticies. 

• Lower  values  of  artificial  dissipation  may  be  used. 

• Larger  time  steps  may  be  used. 

• Improved  convergence  behaviour. 

• More  closely  resembles  experimental  data. 


• Pk>  viewed  as  velocities  in  temporally  parabolized  grid  equa 
tions. 


ADAPTED  DISTRIBUTION  M 
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ADAPTED  GRID 
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GENDE++  - A Multi-Block  Structured  Grid  System 

by 

Tonya  Williams,  Naren  Nadenthiran,  Hugh  Thornburg,  and  Bharat  K.  Soni 
NSF  Engineering  Research  Center  for 
Computational  Field  Simulation 
Mississippi  State  University 
Mississippi  State,  MS  39762 

ABSTRACT 

The  computer  code  GENIE++  (Soni  et  al.  1992)  is  a continuously  evolving  grid  system  containing 
a multitude  of  proven  geometry/grid  techniques.  The  generation  process  in  GEN1E++  is  based  on 
an  earlier  version.  The  process  uses  several  techniques  either  separately  or  in  combination  to  quickly 
and  economically  generate  sculptured  geometry  descriptions  and  grids  for  arbitrary  geometries.  The 
computational  mesh  is  formed  by  using  an  appropriate  algebraic  method.  Grid  clustering  is  accom- 
plished with  either  exponential  or  hyperbolic  tangent  routines  which  allow  the  user  to  specify  a de- 
sired point  distribution.  Grid  smoothing  can  be  accomplished  by  using  an  elliptic  solver  with  proper 
forcing  functions.  B -spline  and  Non-Uniform  Rational  B-splines  (NURBS)  algorithms  are  used 
for  surface  definition  and  redistribution.  The  built-in  sculptured  geometry  definition  with  desired 
distribution  of  points,  automatic  Bezier  curve/surface  generation  for  interior  boundaries/surfaces, 
and  surface  re-distribution  is  based  on  NURBS.  Weighted  Lagrange/Hermite  transfinite  interpola- 
tion methods,  interactive  geometry/grid  manipulation  modules,  and  on-line  graphical  visualization 
of  the  generation  process  are  salient  features  of  this  system,  which  result  in  a significant  time  savings 
for  a given  geometry/grid  application. 

The  development  of  the  system,  as  well  as  computational  examples  of  practical  interest  will  be  pres- 
ented to  demonstrate  the  success  of  these  methodologies.  Complete  documentation  is  available  us- 
ing Mosaic.  Versions  are  available  for  PC’s,  X window,  and  SGI  systems.  It  is  planned  to  place  this 
code  in  the  public  domain  by  August  1995. 
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INITIALIZATION  OPTIONS 

1 TOGGLE  REAL  TIME  PLOTTING 

2 TOGGLE  PROMPTING 

3 GIVE  TITLE  TO  GRID 

4 CHANGE  CURRENT  GRID  BLOCK  SIZE 

5 CHANGE  MAXIMUM  GRID  SIZES 

6 CHANGE  CURRENT  BLOCK  NUMBER 

7 CHANGE  MAXIMUM  NUMBER  OF  BLOCKS 

8 TOGGLE  GRID  GENERATION  MODE 

9 INITIALIZE  DATABASE 

10  INITIALIZE  ZONAL  INFORMATION 

11  VIEW  NON-BLOCK  GRID 

12  VIEU  ONE  BLOCK 

13  VIEW  ALL  BLOCKS 

14  EXIT  INITIALIZATION 

15  QUIT  GRID  GENERATION 
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Surface  and  Volume  grid  generation 
in  parametric  form 

by 

Tzuyi  Yu,  Bharat  K.  Soni 
NSF/ERC  For  Computational  Field  Simulation 
Mississippi  State  University,  MS39762 
Ted  Benjamin,  Robert  Williams 
Marshall  Space  Flight  Center,  AL  35812 


ABSTRACT 


The  algorithms  for  surface  modeling  and  volume  grid  generation  using  parametric  NURBS  geometric 
representation  are  presented.  The  enhanced  re-parameterization  algorithm  which  can  yield  a desired 
physical  distribution  on  the  curve,  surface  and  volume  is  also  presented.  This  approach  bridges  the  gap 
between  CAD  surface  / volume  definition  and  surface  / volume  grid  generation. 

INTRODUCTION 


Surface  grid  generation  is  the  most  labor  intensive  part  of  the  overall  complex  three  dimensional 
grid  generation  process.  Also,  a significant  amount  of  effort  is  required  in  changing  the  resolutions  (grid 
sizes)  and  / or  the  distribution  of  the  grid  while  maintaining  geometry  fidelity.  In  the  last  few  years,  vari- 
ous researchers  have  concentrated  on  utilizing  the  Computer  Aided  Geometry  Design  (CAGD)  techniques 
to  expedite  the  overall  surface  generation  process.  In  this  presentation,  a parametric  formula  which  has 
been  used  in  CAD  system  is  extended  with  re-parameterization  approach  to  numerical  grid  generation  for 
modeling  the  surface  as  well  as  the  volume  grid. 

There  are  many  parametric  approaches  for  representing  sculptured  geometry,  such  as  rational  or 
non-rational  Bezier,  cubic  splines,  rational  or  non-rational  B-spines, ...  ,etc.  Among  these  representation, 
the  Non-Uniform  Rational  B-Splines  (NURBS)  has  been  widely  accepted  among  these  researchers. 
NURBS  has  been  widely  utilized  to  represent  and  design  geometry  in  the  CAD/CAM  and  the  graphics 
community  due  to  its  powerful  features,  such  as  the  local  control  property,  variation  diminishing,  convex 
hull  and  affine  invariance  [Ref  1,2).  Also  the  geometry  tool  kits,  such  as  curve/surface  interpolation,  data 
reduction,  degree  elevation,  knot  insertion  and  splitting,  are  well-developed  [Ref  1,2,].  These  properties 
have  made  NURBS  representation  very  popular  in  recent  developments  in  CAD/CAM.  However,  the  dis- 
tribution requirements  in  CFD  application  are  much  complicated  than  those  in  CAD  system.  Hence,  the 
NURBS  must  be  cooperated  with  re-parameterization  algorithm  so  that  it  can  be  more  useful  in  grid  gen- 
eration. Computational  examples  associated  with  practical  configurations  are  shown  in  Figure  1 and  2. 
The  re-parameterization  approach  described  in  many  research  [Ref  2,3]  is  implemented  by  iteration  pro- 
cess, which  needs  a lot  of  computation  time.  The  more  efficient  and  robust  approach  presented  here  needs 
only  one  interpolation  process. 

The  development  of  the  software  based  on  NURBS  representation  package:  CAGI  (Computer 
Aided  Grid  Interface)  was  initiated  by  authors  under  the  sponsorship  of  NASA  Marshall  Space  Flight 
Center.  The  purpose  of  this  presentation  is  to  present  the  progress  realized  in  enhancing  the  NURBS 
based  curve  / surface  grid  generation  techniques  into  a 3D  volume  grid  generation  technique.  To  this  end, 
various  options  for  generating  3D  volume  geometry-grid  are  discussed.  A reparameterization  scheme 
has  been  developed  to  achieve  desired  distribution  in  physical  space.  Computational  examples  for  model- 
ing practical  configurations  have  been  exercised  using  the  volume  options  and  the  reparameterization 
scheme. 
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Figure  1 : Propulsional  example. 

3D  NURBS  control  patches  model 
the  missile  (with  fins)  geometry. 


Figure  2: 

3D  NURBS  control  patches  model 
the  single  rotation  propfan. 
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CAGI : A NURBS  database  package 
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NURBS  control  volume  for  a nozzle  with  square-ellipse-circle  sections 

( O Type  Control  Volume ) 
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Overview  of  CFD  Analyses  Supporting  the  Reusable  Solid  Rocket 

Motor  (RSRM)  Program  at  MSFC 

E.  Stewart,  P.  McConnaughey,  J.  Lin,  E.  Reske,  and  D.  Doran,  NASA/MSFC,  R.  H. 

Whitesides,  ERC,  Inc.,  Huntsville,  AL,  and  Y.-S.  Chen  ESI,  Huntsville,  AL 

During  the  past  year,  various  CFD  analyses  were  performed  at  MSFC  to 
support  the  RSRM  program.  The  successful  completion  of  these  analyses 
was  realized  through  the  cooperation~df  ESI,  ERC,  and  The  Computational 
Fluid  Dynamics  Branch  (ED32)  at  MSFO~and  involved  application  of  the  CFD 
codes  FDNS  and  CELMINT.  The  topics  addressed  by  the  analyses  were;  1. 
the  design  and  prediction  of  slag  accumulation  within  the  five  inch  test 
motor,  2.  prediction  of  slag  pool  behavior  and  its  response  to  lateral 
accelerations,  3.  the  clogging  of  potential  insulation  debonds  within  the 
nozzle  by  slag  accumulation,  4.  the  behavior  of  jets  within  small  voids 
inside  nozzle  joint  gaps,  5.  the  effect  of  increased  inhibitor  stiffness  on 
motor  acoustics,  and  6.  the  effect  of  a nozzle  defect  on  particle 

impingement  enhanced  erosion.  Topics  1,  2,  and  5 will  be discussed  in 

some  detail  by  other  speakers  at  the  conference  and  are  only  mentioned 
here  -for  the  sake  of  completeiie^srThus^lie  emphasis  of  this  presentation 
will  be  to  further  discuss  the  work  involved  in  topics  3,  4,  and  6. 

Eric  Stewart 
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Insulation  debond  analysis 

Potential  RTV  flaw  analysis 

Nose  inlet  assembly  wetline  Investigation 
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impact  m joint  gap  clogging  prediction  methodology  is  available  to 

support  potential  anomalies 
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• Results  (Problem  1,  debond  vent  to  ambient) 


gap  width 

0.010" 

0.005" 

0.002" 


time  to  clog 

0.05  sec. 
0.02  sec. 
0.006  sec. 


• Results  (Problem  2,  start-up  transient  through  debond) 

- clogging  of  debond  predicted  in  0.61  sec.  after  initiation  of  particle 
flow 

- lower  mass  flow  rate  (4X  less)  due  to  cavity  fill  results  in  fewer 
particles  to  condense  on  pore  wall 
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Joint  Gap  Clogging  Prediction 
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(rather  than  1.25") 

Jet  spreading  predictions  are  available  to  support  potential 
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Nozzle,  Joint  #4  Computational  Geometry 
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George  C.  Marshall  Space  Flight  Cer 


Velocity  Magnitudes 


Center  Line  Velocity  Profile 
Through  Comer  A 


Average  Velocity  Magnitudes 
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Nose  Inlet  Assembly  Wetline  Investigation 


Issues  • Enhancement  of  nozzle  erosion  due  to  presence  of  defect 

• Effect  of  defect  on  slag  particle  impingement 

Approach  • Assume  wedge  shaped  nozzle  defect 

• Use  two-phase  flow  results  to  assess  flow  environment 
near  defect 

• Use  current  data/ experience  base  to  assess  potential  flow 
deviations 

Size  of  defect  relative  to  local  boundary  layer  is  not  sufficient 
to  significantly  alter  flow  external  to  boundary  layer 
Main  source  of  particle  impingement  is  external  to  boundary 
layer 

Erosion  enhancement  due  to  particle  impingement  is  not 
significantly  altered  by  presence  of  defect 

Impact  • Recommend  nozzle  in  question  for  flight 


Results  • 
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Geometry  of  Nozzle  Nose  Region 


ED32  Results 
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Improve  readiness  to  address  potential  anomalies 
- Perform  similar  analyses  at  additional  burn  times 


RSRM  Chamber  Pressure  Oscillations: 
Transit  Time  Models  and  Unsteady  CFD 


Tom  Nesman  and  Eric  Stewart 
Fluid  Dynamics  Division 

National  Aeronautics  and  Space  Administration,  MSFC 
Marshall  Space  Flight  Center,  AL  35812 


Abstract 

Space  Shuttle  solid  rocket  motor  (SRM)  low  frequency  internal  pressure  oscillations  have  been 
observed  since  early  testing.  The  same  type  of  oscillations  are  also  present  in  the  redesigned  solid 
rocket  motor  (RSRM).  The  oscillations,  which  occur  during  RSRM  bum,  are  predominantly  at  the 
first  three  motor  cavity  longitudinal  acoustic  mode  frequencies.  Broadband  flow  and  combustion 
noise  provide  the  energy  to  excite  these  modes  at  low  levels  throughout  motor  burn,  however,  at 
certain  times  during  bum  the  fluctuating  pressure  amplitude  increases  significantly.  The  increased 
fluctuations  at  these  times  suggests  an  additional  excitation  mechanism. 

The  RSRM  has  inhibitors  on  the  propellant  forward  facing  surface  of  each  motor  segment.  The 
inhibitors  are  in  a slot  at  the  segment  field  joints  to  prevent  burning  at  that  surface.  The  aft  facing 
segment  surface  at  a field  joint  slot  bums  and  forms  a cavity  of  time  varying  size.  Initially  the 
inhibitor  is  recessed  in  the  field  joint  cavity.  As  propellant  burns  away  the  inhibitor  begins  to 
protrude  into  the  bore  flow.  Two  mechanisms  (transit  time  models)  that  are  considered  potential 
pressure  oscillation  excitations  are  cavity  edge-tones,  and  inhibitor  hole-tones.  Estimates  of 
frequency  variation  with  time  of  longitudinal  acoustic  modes,  cavity  edge-tones,  and  hole-tones 
compare  favorably  with  frequencies  measured  during  motor  hot  firing.  It  is  believed  that  the 
highest  oscillation  amplitudes  occur  when  vortex  shedding  frequencies  coincide  with  motor 
longitudinal  acoustic  modes. 

A time  accurate  CFD  analysis  was  made  to  replicate  the  observations  from  motor  firings  and  to 
observe  the  transit  time  mechanisms  in  detail.  FDNS  is  the  flow  solver  used  to  detail  the  time 
varying  aspects  of  the  flow.  The  fluid  is  approximated  as  a single-phase  ideal  gas.  The  CFD 
model  was  an  axisymmetric  representation  of  the  RSRM  at  80  seconds  into  burn.  Deformation  of 
the  inhibitors  by  the  internal  flow  was  determined  through  an  iterative  structural  and  CFD  analysis. 
The  analysis  domain  ended  just  upstream  of  the  nozzle  throat.  This  is  an  acoustic  boundary 
condition  that  caused  the  motor  to  behave  as  an  closed-open  organ  pipe.  This  differs  from  the 
RSRM  which  behaves  like  a closed-closed  organ  pipe. 

The  unsteady  CFD  solution  shows  RSRM  chamber  pressure  oscillations  predominantly  at  the 
longitudinal  acoustic  mode  frequencies  of  a closed-open  organ  pipe.  Vortex  shedding  in  the  joint 
cavities  and  at  the  inhibitors  contribute  disturbances  to  the  flow  at  the  second  longitudinal  acoustic 
mode  frequency.  Further  studies  are  planned  using  an  analysis  domain  that  extends  downstream  of 
the  nozzle  throat. 
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RSRM  - Chamber  Pressure  Oscillations: 

Transit  Time  Models  and  Unsteady  CFD 

Workshop  for  CFD  Applications  in  Rocket  Propulsion  and  Launch 

Vehicle  Technology 


Tom  Nesman  and  Eric  Stewart 
Fluid  Dynamics  Division  - NASA  - MSFC 


RSRM  Pc  Oscillations;  Transit  Time  George  C.  Marshall  Space  Flight  Center 

Models  and  Unsteady  CFD  Fluid  Dynamics  Division 
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Workshop  for  CFD  Applications  in  Rocket  Propulsion  and  Launch  Vehicle  Technology 


RSRM  Pc  Oscillations  George  C.  Marshal!  Space  Flight  Center 

Fluid  Dynamics  Division 
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Typical  RSRM  Pc  Isoplot 


George  C.  Marshall  Space  Flight  Center 
Fluid  Dynamics  Division 


Head-end  Chamber  Pressure  (Pc)  Measurement 
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Excitation  Mechanisms 


Cavity  Edge-Tone  Mechanism 


Acoustic  waves-" 

- 

I;!;,-' 


Downstream 

impingement 


Feedback  Transit  Time  Models 

Empirical  estimates  - time  for  disturbance  to 
travel  downstream  plus  time  for  feedback  to 
sensitive  region 

Analytical  estimates  - sum  of  phase  of 
downstream  traveling  wave  and  phase  of 
upstream  traveling  acoustic  wave 


\ ' Field  joint  cavity 

Oscillating  shear  layer 


Sensitive  region 


Upstream 
inhibitor  * 


Inhibitor  to  Inhibitor  Hole-tone  Mechanism 
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RSRM  Unsteady  CFD  Timeplot 


Fluctuating  Pressure  Downstream  of  Inhibitor  Tip 
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RSRM  Pc  Oscillations:  Transit  Time  George  c.  Marshall  Space  Flight  Center 

Models  and  Unsteady  CFD  Fluid  Dynamics  Division 
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A Coupled  CFD/FEM  Structural  Analysis  to  Determine  Deformed  Shapes  of 

the  RSRM  Inhibitors 


Richard  A.  Dill  and  R.  Harold  Whitesides 
ERC,  Incorporated,  Hunstville,  AL  35816 


Recent  trends  towards  an  increase  in  the  stiffness  of  the  NBR  insulation  material  used  in 
the  construction  of  RSRM  propellant  inhibitors  prompted  questions  about  possible  effects 
on  RSRM  performance.  The  specific  objectives  of  the  CFD  task  included:  1)  the 
definition  of  pressure  loads  to  calculate  the  deformed  shape  of  stiffer  inhibitors,  2)  the 
calculation  of  higher  port  velocities  over  the  inhibitors  to  determine  shifts  in  the  vortex 
shedding  or  edge  tone  frequencies  and  3)  the  quantification  of  higher  slag  impingement 
and  collection  rates  on  the  inhibitors  and  in  the  submerged  nose  nozzle  cavity. 

A coupled  CFD/Finite  element  structural  analysis  was  required  to  calculate  the  deformed 
inhibitor  geometry.  Since  the  NBR  inhibitor  material  erodes  at  a different  rate  than  the 
motor  propellant  burns,  an  inhibitor  stub  which  protrudes  above  the  propellant  into  the 
port  cavity  is  created  during  motor  operation.  The  impinging  port  flow  causes  the 
inhibitor  stub  to  bend  in  the  downstream  flow  direction.  Since  a stiffer  NBR  inhibitor 
material  would  cause  the  inhibitor  to  bend  less,  it  was  necessary  to  know  the  difference 
in  the  bending  of  the  original  NBR  material  compared  to  the  stiffer  NBR  materiai.  The 
CELMINT  CFD  computer  code  was  used  to  perform  the  fluid  dynamic  calculations  of  the 
motor  flow  field.  The  structural  bending  effect  of  the  pressure  loads  from  the  CFD  code 
was  analyzed  by  ED28.  Initially,  the  CELMINT  code  was  used  to  determine  the  flow  field 
and  inhibitor  pressure  loads  for  unbent  motor  inhibitors.  This  pressure  loading  on  the 
inhibitors  was  used  by  ED28  to  generate  the  bending  which  would  occur  in  the  inhibitor. 
The  computed  bent  inhibitor  geometry  was  then  used  again  by  the  CFD  code  to  compute 
a new  pressure  loading  on  the  inhibitors.  This  iterative  computation  between  the  CFD 
code  and  the  structural  analysis  code  was  continued  until  convergence  in  the  inhibitor 
bent  geometry  was  achieved. 

The  CFD  solution  was  then  used  to  assess  the  effect  of  higher  flow  velocities  and  edge 
tone  frequencies  from  the  reduced  inhibitor  bending  on  the  maximum  oscillating  pressure 
amplitudes  that  occur  during  resonance  between  the  edge  tones  and  the  motor 
longitudinal  modes.  Also,  a comparison  of  the  difference  in  slag  accumulation  between 
the  two  NBR  materials  was  also  made  to  determine  if  the  stiffer  material  increases  slag 
collection  in  the  field  joints  and  the  submerged  nozzle  cavity. 

The  coupled  CFD/FEM  structural  analysis  was  successful  in  defining  the  effect  of 
inhibitor  stiffness  on  inhibitor  geometry  and  the  shift  in  edge  tone  frequencies.  Also,  the 
two-phase  CFD  analysis  showed  that  there  was  a small  increase  in  the  rate  of  slag 
accumulation  at  the  aft  inhibitor;  however,  motor  trajectory  analyses  of  slag  debris  shed 
from  the  inhibitors  showed  that  the  debris  would  pass  out  the  motor  nozzle  and  therefore 
create  no  additional  slag  accumulation  in  the  slag  pool  around  the  nozzle. 
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Background 


• In  October,  1994,  Thiokol  reported  the  use  of  NBR  material  in  RSRM’s  with  properties 
significantly  different  from  the  historical  database. 

• A 30%  to  40%  increase  in  modulus  was  reported. 

• This  increased  stiffness  had  the  potential  to  affect  the  amplitude  of  chamber  pressure 
oscillations  in  the  SRM: 

-By  changing  the  inhibitor  structural  response 

-By  indirectly  changing  the  flow/acoustic  interaction 

• The  slag  accumulation  in  the  field  joints  and  submerged  nozzle  region  might  also  be 
increased  thereby  increasing  the  potential  for  pressure  and  thrust  perturbations. 


ERC,  Inc. 
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Objectives  of  Coupled  CFD/FEM  and  Two-Pbase  CFD  Analyses 


• Determine  deformed  geometry  of  NBR  inhibitors  at  the  forward,  center  and  aft  joints  for 
both  nominal  and  stiff  NBR  materials  using  a coupled  CFD/FEM  analysis. 

® Determine  effect  of  inhibitor  properties/geometry  on  inhibitor  hole  velocities  to  evaluate 
effect  on  hole  edge  tone  (vortex  shedding)  frequencies. 

© Determine  effect  of  inhibitor  properties/geometry  on  slag  accumulation  on  both  the 
inhibitor  surfaces  and  underneath  the  nozzle  nose. 


ERC,  Inc. 


Coupled  CFD/FEIVI  Analysis  Approach 


1 ) Perform  single-phase  gas  CFD  analysis  of  entire  RSRM  port  at  80  second  burn  time 
using  straight  inhibitor  lengths  from  erosion  analysis. 

2)  Perform  FEM  structural  analysis  on  inhibitors  to  determine  deformations  using  surface 
pressure  distributions  from  CFD  analysis. 

3)  Perform  CFD  analysis  using  deformed  inhibitor  geometries  from  step  2). 

4)  Repeat  steps  2)  and  3)  until  convergence  of  inhibitor  geometry  is  achieved. 

5)  Provide  velocity  profile  at  each  inhibitor  location  for  both  nominal  and  stiff  inhibitors  as 
input  to  flow/acoustic  interaction  analysis. 


ERC,  Inc. 


Two-Phase  Flow  CFD  Methodology 
CELMINT  Code 

(Combined  Euferian  Lagrangian  Multi-Dimensional  Implicit  Nonlinear  Time-Dependent) 


• Navier-Stokes  Solution 

- Fully  implicit,  density-based,  conservative,  ensemble-averaged  Navier-Stokes  code 

- Low  and  high  Reynolds  number  and  wall  injection  k-e  models 

- Equilibrium  and  finite-rate  chemistry  for  multi-species  flows 


• Two-phase  Flow  Models 

- Coupled  Eulerian-Lagrangian  for  solid  and  liquid  phases 

- Hermsen  aluminum  burn  rate  model  for  particle  combustion 

- Specification  of  particle  properties  (density,  size  distribution) 

- Particle  break-up  based  on  Weber  number 

- Agglomeration  based  on  collisions  between  discrete  phase  particles  and  continuous 
phase  smoke  particles 

- Programmable  for  various  particle  capture  criteria 


ERC,  Inc. 


Propellant  Thermochemical  Properties  and  Motor  Operating  Conditions 

RSRM  80  Second  Burn  Time 
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RSRM  Motor  Geometry 
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(Dashed  Line  Shows  The  67  Second  Bum  Back) 


Computational  Grid  Resolution 
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RSRM  80  Second  Stiff  NBR  Inhibitor 
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RSRM  80  Second  Stiff  NBR  Inhibitor 
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Computational  Grid,  Aft  Slot 

RSRM  80  Second  Stiff  NBR  Inhibitor 
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RSRM  80  Second  Stiff  NBR  Inhibitor 
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Velocity  Vectors,  Aft  Slot 

RSRM  80  Second  Stiff  NBR  Inhibitor 


Velocity  Vectors,  Submerged  Nozzle 

RSRM  80  Second  Stiff  NBR  Inhibitor 


Forward  Inhibitor  Radial  Pressure  Distribution 
RSRHI  80  Seconds  Burn  Time 
Stiff  HER 
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iiotor  Radius  finches) 


Center  Inhibitor  Radial  Pressure  Distribution 
RSRM  §0  Seconds  Burn  Time 
Stiff  NBR 
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Motor  Radius  (inches) 


Aft  inhibitor  Radial  Pressure  Distribution 
RSRM  80  Seconds  Burn  Time 
Stiff  NBR 
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Motor  Radius  finches) 


Forward  Inhibitor  Deformation  iterations 

Stiff  HER 

RSRM  80  Seconds  Burn  Time 
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852  856  860  864  868  872 

Axial  Motor  Station  (Inches) 


Forward  Inhibitor  Deformations 
Nominal  and  Stiff  NBR 
RSRM  80  Seconds  Burn  Time 
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Axial  Motor  Station  (inches) 


Center  Inhibitor  Deformations 
Nominal  and  Stiff  NBR 
RSRIi  80  Seconds  Burn  Time 
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Axial  Motor  Station  (inches) 


Aft  Inhibitor  Deformations 
Nominal  and  Stiff  NBR 
RSRM  80  Seconds  Burn  Time 
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Axial  Motor  Station  finches) 


Port  Velocity  Profile  at  Forward  Inhibitor 
Nominal  and  Stiff  NBR 
RSRsi  80  Seconds  Burn  Time 
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Axial  Velocity  (ft/sec) 


Port  Velocity  Profile  at  Center  Inhibitor 
Nominal  and  Stiff  HER 
RSRfU  80  Seconds  Burn  Time 
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Axial  Velocity  (fft/sec) 


Port  Velocity  Profile  at  Aft  Inhibitor 
Nominal  and  Stiff  NBft 
RSRM  80  Seconds  Burn  Time 
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100  150  200  250  300  350  400  450  500  550 

Axial  Velocity  (ft/sec) 


Comparison  of  the  Motor  Port  Velocity  Profiles  immediately 

Upstream  ©f  Wozzle  Wose 
RSRii  80  Seconds  Burn  Time 
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Coupled  CFD/FEM  Analysis  Conclusions 
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Two-Phase  CFD  Analysis  Approach 


• Perform  two-phase  CFD  analysis  of  RSRM  port  at  80  second  burn  time  using  final 
deformed  inhibitor  geometries  for  both  nominal  and  stiff  inhibitors. 

• Calculate  slag  captured  on  both  nominal  and  stiff  inhibitors  at  all  three  field  joints. 

• Perform  trajectory  analysis  for  slag  debris  shed  from  inhibitor  tips  for  all  above  cases  to 
determine  whether  it  passes  through  nozzle  or  accumulates  underneath  nozzle  nose. 


ERG,  Inc . 


RSRM  Inhibitor  Slag  Accumulation 
Nominal  and  Stiff  NBR  Inhibitors 
SCI  Second  Burn  Time 
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ectory  Results 
NBR  Inhibitors 
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Two=Phas@  CFD  Analysis  Conclusions 


• The  rate  of  slag  accumulation  for  both  the  nominal  and  stiff  inhibitors  at  all  joints  is  a very 
small  percentage  of  the  total  motor  slag  accumulation  rate. 

• The  rate  of  slag  accumulation  on  the  center  inhibitor  is  approximately  four  times  greater 
for  the  stiff  NBR  compared  to  the  nominal  NBR. 

• Slag  debris  shed  from  the  nominal  inhibitors  at  all  three  joints  exits  the  nozzle  throat 
plane. 

• Slag  debris  shed  from  the  stiff  inhibitors  at  the  forward  and  center  joints  exits  the  nozzle 
throat  plane.  Slag  from  the  aft  joint  stiff  inhibitor  impacts  the  nozzle  entrance  ramp. 

• No  excess  slag  collected  on  the  stiff  inhibitors  is  transported  underneath  the  nozzle  nose 
to  add  to  the  normal  slag  pool. 
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NASA,  Marshall  Space  Flight  Center 
Computational  Fluid  Dynamics  Branch 
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Summary  and  Future  Plans 
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--  appropriate  methodology? 
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4 AXIAL  STATIONS 


ASRM  Air  Fiow  Model  Inlet  Velocity  Profiles  - Redesigned  Velocity 
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CONTOUR  Itt/ELS 
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Application  of  Two-Phase  CFD  to  the  Design  and  Analysis  of  a Subscale 
Motor  Experiment  to  Evaluate  Propellant  Slag  Production 

R.  Harold  Whitesides  and  Richard  A.  Dili 
ERC,  Incorporated,  Huntsville,  AL  35816 


Abstract 

The  RSRM  Pressure  Perturbation  Investigation  Team  concluded  that  the  cause  of  recent 
pressure  spikes  during  both  static  and  flight  motor  burns  was  the  expulsion  of  molten 
aluminum  oxide  slag  from  a pool  which  collects  in  the  aft  end  of  the  motor  around  the 
submerged  nozzle  nose  during  the  last  half  of  motor  operation.  It  is  suspected  that  some 
motors  produce  more  slag  than  others  due  to  differences  in  aluminum  oxide  agglomerate 
particle  sizes  which  may  relate  to  subtle  differences  in  propellant  ingredient 
characteristics  such  as  particle  size  distributions,  contaminants,  or  processing  variations. 

In  order  to  determine  the  effect  of  suspect  propellant  ingredient  characteristics  on  the 
propensity  for  slag  production  in  a real  motor  environment,  a subscale  motor  experiment 
was  designed  to  accomplish  this  objective.  An  existing  5 inch  ballistic  test  motor  was 
selected  as  the  basic  test  vehicle  due  to  low  cost  and  quick  turn  around  times.  The 
standard  converging/diverging  nozzle  was  replaced  with  a submerged  nozzle  nose 
design  to  provide  a positive  trap  for  the  slag  which  would  increase  both  the  quantity  and 
repeatability  of  measured  slag  weights.  CFD  was  used  to  assess  a variety  of  submerged 
nose  configurations  to  identify  the  design  which  possessed  the  best  capability  to  reliably 
collect  slag.  CFD  was  also  used  to  assure  that  the  final  selected  nozzle  design  would 
result  in  flow  field  characteristics  such  as  dividing  streamline  location,  nose  attach  point, 
and  separated  flow  structure  which  would  have  similitude  with  the  RSRM  submerged 
nozzle  nose  flow  field.  It  was  also  decided  to  spin  the  5 inch  motor  about  its  longitudinal 
axis  to  further  enhance  slag  collection  quantities.  Again,  CFD  was  used  to  select  an 
appropriate  spin  rate  along  with  other  considerations,  including  the  avoidance  of  burn 
rate  enhancement  from  radial  acceleration  effects. 


The  CFD  analyses  were  performed  with  the  CELMINT  code  which  is  a two-phase 
Navier-Stokes  coded  employing  an  Eulerian/Lagrangian  scheme,  a low  Reynolds 
number  k-s  turbulence  model  modified  for  wall  injection,  and  both  surface  and  distributed 
particle  combustion  models  which  include  particle  agglomeration  and  break-up. 
Aluminum  oxide  particle  distributions  were  measured  with  RSRM  propellant  in  a 
combustion  bomb  with  particle  quench  capability.  Predictions  for  slag  weights  and  slag 
distribution  patterns  were  compared  with  slag  weight  data  from  defined  zones  in  the 
motor  and  nozzle.  Various  parameters  were  investigated  to  reconcile  differences 
between  CFD  predictions  and  data.  General  comparisons  were  acceptable  considering 
combustion  bomb  data  on  particle  sizes  was  not  available  for  each  propellant  sample. 
Confidence  in  using  this  methodology  in  the  RSRM  was  enhanced  by  this  successful 
subscale  experiment. 
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Background 


Flight  and  static  test  data  for  the  Space  Shuttle  Reusable  Solid  Rocket  Motors  reveals 
roughness  and  small  spikes  in  the  pressure  trace  for  some  motors  during  the  65-75 
second  time  period. 
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Experimental  Program  Objective  and  Approach 


Objective: 

Develop  and  employ  a subscale  rocket  test  motor  capable  of  measuring  relative  slag 

production  of  propellants  with  subtle  changes  in  ingredient  characteristics. 

Approach; 

• Use  an  existing  Thiokol  5-inch  diameter  ballistic  test  motor  and  static  test  spin  stand. 

• Modify  existing  converging/diverging  nozzle  entrance  geometry  by  incorporating 
submerged  nose  to  enhance  slag  capture  and  retention. 

• Select  motor  operating  pressure  to  match  full  scale  motor  pressure.  Select  spin  rate  to 
enhance  slag  capture  but  avoid  propellant  burn  rate  augmentation. 

• Use  CFD  to  determine  overall  viability  of  experiment,  to  aid  in  design  of  motor 
components,  to  support  selection  of  test  conditions,  and  to  analyze  test  results. 
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Specific  CFD  Analysis  Tasks 
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Use  analysis  results  to  upgrade  two-phase  CFD  model  for  RSRM  slag  predictions. 


SPIN  MOTOR 
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Head  End  Closure 


FIVE  INCH  SPIN  MOTOR  NOZZLE 

SLAG  ZONES 
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PIN  MOTOR 
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Velocity  Vectors,  Submerged  Nozzle 
RSRM  80  Second  Stiff  NBR  Inhibitor 


Velocity  Vectors,  Submerged  Nozzle 
Spin  Motor  at  50%  Web  Tim© 


Two-Phase  Flow  CFD  Methodology 
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Programmable  for  various  particle  capture  criteria 


Effect  of  Spin  Rate  on  Nozzle  Slag  Accumulation 
WECCO  AP,  Surface  Combustion,  400RPM 
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Spin  Rate  (RPM) 


Computational  Grid  Resolution 
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Spin 
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Spin  Motor  At  50  % Web  Time 
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Spin  Motor  At  85  % Web  Time 


RSRM  Propellant  Thermochemical  and  Nominal  Particle  Properties 
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Spin  Motor  Operation  Conditions 


Parameter 
Chamber  Pressure: 
Mass  Flow  Rate: 
Propellant  Burning  Area: 

Throat  Diameter: 

Burn  Time: 


15%  Web 
610.6  psia 
2.613  Ibm/s 
115.47  in2 


50%  Web 
628.8  psia 
2.691  Ibm/s 
114.67  in2 


85%  Web 
610.6  psia 
2.613  Ibm/s 
112.35  in2 


.916  inches 
2.71  seconds 

ERC,  Inc. 


Experimental  Data  and  Curve  Fit  Distribution  Functions 
WECCO  Quench  Bomb  Data™  500  psi 
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Particle  Diameter  {microns) 


Experimental  Data  and  Curve  Fit  Density  Functions 
WECCO  Quench  Bomb  Data-  500  psi 
3-inch  Quench  Distance 


uojpejj  sseyv 
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Particle  Diameter  (microns) 


Experimental  Data  and  Curve  Fit  Distribution  Functions 
Kerr  McGee  Quench  Bomb  Data-  500  psi 
3-inch  Quench  Distance 


n n 
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Particle  Diameter  (micron®) 


Experimental  Data  and  Cyrve  Fit  Density  Functions 

Kerr  iWlcGe©  Quench  Bomb  Data-  500  psi 
3-inch  Quench  Distance 


uojpejj  sseyy 
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Particle  Diameter  (microns) 


Effect  of  Particle  Size  Distribution  on  Slag  Accumulation 
WECCO  AP,  Surface  Combustion,  400  rpm 


1267 


Percent  RHotor  Web  Time 


Nozzle  Slag  Accumulation 
WECCO  and  Kerr  fJcGee  AP 
Surface  Combustion,  400  RPiS 
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Percent  motor  Web  Time 
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Spin  Motor  at  50%  Web  Time 
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Velocity  Field  In  The  Submerged  Nozzle  Region 
Spin  Motor  At  15  % Web  Time 
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pin  Motor  At  50  % Web  Time 
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Velocity  Field  In  The  Submerged  Nozzle  Region 

Spin  Motor  At  85  % Web  Time 
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Spin  Motor  At  85  % Web  Time 
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Spin  Motor  At  85  % Web  Time 
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Spin  liotor  At  85  % Web  Time 


Spin  Motor  Circumferential  Velocity 
WECCO,  Surface  Combustion,  400  RPM 

15%  Web  Time 
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Circumferential  Velocity  (feet/second) 
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Nozzle  Slag  Accumulation  Rate 
WECCQ  APS  Surface  Combustion,  400  rpm 
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Percent  Motor  Web  Time 


Nozzle  Slag  Accumulation 
WECCO  AP,  Surface  Combustion,  400  rptn 
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Percent  Motor  Web  1 


Slag  Accumulation  Per  Unit  Area  Along  The  Nozzle 
WECCO  AP,  Surface  Combustion,  400  rpm 


ngth  Along  Nozzle  Surface  From  The  Throat  fin. 
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of  Initial  Particle  Velocity  on  Slag  Accumulation 
ftfECCO  AP,  Surface  Combustion,  400  rpm 
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Percent  Motor  Web  Time 


Effect  of  Particle  Density  on  Slag  Impingement  Rates 
WECCO  AP,  Surface  Combustion,,  400  RPli 

85%  Web  Time 
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Nozzle  SSag  Zones 


Nozzle  Slag  Accumulation 
WECCO  AP,  Distributed  Combustion,  400  rpm 
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Percent  Motor  Web  Time 
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Numerical  Investigation  of  Slag  Behavior  for  R8RM 

P.  Liaw  ,Y-S  Chen,  H.  Shang,  amd  M.  Shih 
Engineering  Sciences,  Inc.,  Huntsville,  AL 
and 

D.  Doran  and  E.  Stewart 
NASA  MSFC,  Huntsville,  AL 

ABSTRACT 

A ‘>  A/a ^ ’ 

It  is  known  that  the  flowfield  of  the  SRM  (Solid  Rocket  Motor)  is  very  complicated  due  to  the 
complex  characteristics  of  turbulent  multi-phase  flow,  chemical  reaction,  particle  combustion,  evaporation, 
breakup  and  agglomeration  etc.  It  requires  multi-phase  calculations,  chemical  reaction  simulation,  and  particle 
combustion,  evaporation,  and  breakup  models  to  obtain  a better  understanding  of  thermophysics  for  the  SRM 
design  using  numerical  methods.  Also,  the  slag  buildup  due  to  the  molten  particles  is  another  factor  affecting 
the  performance  of  the  SRM.  Thus,  a more  realistic  simulation  is  needed  to  provide  a better  design  guide  to 
improve  the  performance  of  SRM.  To  achieve  this  goal,  the  VOF  (Volume  Of  Fluid)  method  is  used  to 
capture  the  free  surface  motion  so  as  to  simulate  the  accumulation  of  the  molten  particles  (slag)  of  SRM.  A 
Finite-rate  chemistry  model  is  used  to  simulate  the  chemical  reaction  effects.  For  multi-phase  calculations, 
Hermsen  combustion  model  is  used  for  the  AL  particle  combustion  analysis  and  Taylor  Analogy  Breakup 
(TAB)  model  is  used  for  the  particle  breakup  analysis.  An  interphase  mass-exchange  model  introduced  by 
Spalding  is  used  for  the  evaporation  calculation.  The  particle  trajectories  are  calculated  using  a one-step 
implicit  method  for  several  groups  of  particle  sizes  by  which  the  drag  forces  and  heat  fluxes  are  then  coupled 
with  the  gas  phase  equations. 

The  preliminary  results  predicted  a reasonable  physical  simulation  of  the  particle  effects  using  a 
simple  2-D  solid  rocket  motor  configuration.  It  shows  that  the  AL/AL203  particle  sizes  are  reduced  due  to 
the  combustion,  evaporation,  and  breakup.  The  flowfield  is  disturbed  by  the  particles.  Mach  number 
distributions  in  the  nozzle  are  deformed  due  to  the  effect  of  particle  concentrations  away  from  the  center  line. 

The  RSRM  (Redesigned  Solid  Rocket  Motor)  geometry  at  67  seconds  is  employed  to  investigate  the 
slag  behavior  in  the  aft-end  cavity  with  the  combustion,  evaporation,  and  breakup  models.  The  particulate 
phase  was  assumed  to  be  aluminum  oxide  (AL203)  for  the  preliminary  study.  It  is  assumed  that  the 
propellant  grain  of  the  aft-end  cavity  has  burned  out  completely  at  67  seconds.  The  geometry  and  mass  flow 
rate  information  were  provided  by  the  NASA  Marshall  Space  Flight  Center.  The  slag  may  flow  out  of  the 
cavity  and  enter  the  nozzle  due  to  the  accelerations.  The  molten  particles  entering  the  aft-end  cavity  merge 
with  the  slag.  The  volume  of  the  slag  will  grow  and  affect  the  performance  of  the  RSRM.  This  shows  that  the 
effects  of  particles  and  slag  on  the  flowfield  are  very  significant.  From  the  caculation,  a flow  vortex  exists  in 
the  aft-end  cavity  of  the  RSRM.  A stagnation  point  on  the  wall  is  captured.  This  flow  impingement  may 
cause  the  erosion  of  the  wall.  The  shape  of  the  vortex  is  changed  due  to  the  slag.  The  particles  entering  the 
cavity  may  become  slag  and  either  flow  into  or  out  of  the  cavity  depending  on  the  temperature  and  the  surface 
tension  of  the  molten  particles.  An  axial  gravity  force  of  2.4g  is  assumed  to  simulate  the  RSRM  flowfield  at 
67  seconds. 

The  flowfield  analysis  using  the  FDNS  code  in  the  present  research  using  the  proposed  models  should 
provide  a design  guide  for  the  solid  rocket  motors.  The  obtained  results  can  give  the  designer  a basic  guide 
line  for  the  use  of  materials  and  the  nozzle  geometry  to  improve  the  performance  of  SRMs. 
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BACKGROUND  & GENERAL  APPROACH 

• ACCUMULATED  SLAG  WAS  FOUND  IN  THE  AFT-END 
CAVITY  AND  NOZZLE. 

==>  WILL  THIS  AFFECT  THE  MOTOR  PERFORMANCE  DUE 
TO  ITS  EFFECT  ON  THE  PRESSURE? 

• WILL  VOF  METHOD  WORK  FOR  THE  ANALYSIS  OF  SLAG 
BEHAVIOR? 


FORWARD  DOME 
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VOF  Model 

The  V OF  transport  equation  is  given  below: 


da 

~a 


+ 


where  a = 1 stands  for  liquid  and  a = 0 is  for  gas. 
The  interface  is  located  at  1 > a > 0.  For  a given 
solution  of  a field,  equation  (6)  can  be  recast  as: 


for  compressible  gas : 


+ 


dp\u  - ug ).  </> 


3x  ■ 


S ^ , a < 0.01 


for  incompressible  gas : 


d</>  ( \ d(j) 

pM  ~^+  pAu~us)l^~=s*’a  -°-° 

and 

Pm  = Max[pg,ap(  j 


The  interface  a solution  compression  procedure  is 
expressed  as: 


new 


Maxi 0,  Min  1 , 0.5  + f[a 


old 
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and 


(. Interface  volume j 


new 


(. Interface  volume ) 


initial 


The  surface  tension  forces  in  the  continuum 
surface  force  model  is  fomiulated  as  continuous 
body  forces  across  the  interface.  These  forces  can 
be  written  as: 


F_ 


a V n 


) 


a 


A 


\ 


f 


F .=  - a V n 


a y + 


a 


v y j 


, for  2D , axisym  m etric  case  only 


r 


Fz  = - <j|^V  n ja  z , for  3 D case  only 


w here 


surface  tension  constant 


A A 


Vn  = a xx  + ccyy  + a 


a is  0.5  for  the  free  surface.  The  VOF  method  is 
used  to  represent  the  tracking  of  the  free  surface 
between  the  liquid  and  gas  phase. 
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APPLICATIONS 
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3 Wm 


Mach  number  contours  of  RSRM,  no  particles. 
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Velocity  vectors  near  the  aft-end  cavity,  no  particles. 
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Particle  trajectories. 
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Slag  Accumulation  for  3D  RSRM 
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Accumulation  for  3D  RSRM  at 
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Slag  Accumulation  for  3D  RSRM  at  T = 0.0316 
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lag  Accumulation  for  3D  RSRM  at 
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Combustion  Processes  in  Hybrid  Rocket  Engines 

S.  Venkateswaran  and  C.  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802. 


^ p 

In  recent  years,  there  has  been  a resurgence  of  interest  in  the  development  of 
hybrid  rocket  engines  for  advanced  launch  vehicle  applications.  Hybrid  propulsion 
systems  use  a solid  fuel  such  as  hydroxyl-terminated  polybutadiene  (HTPB)  along 
with  a gaseous/liquid  oxidizer.  The  perfomance  of  hybrid  combustors  depend  on  the 
convective  and  radiative  heat  fluxes  to  the  fuel  surface,  the  rate  of  pyrolysis  in  the  solid- 
phase  and  the  turbulent  combustion  processes  in  the  gaseous-phase.  These  processes  in 
combination  specify  the  regression  rates  of  the  fuel  surface  and  thereby  the  utilization 
efficiency  of  the  fuel.  In  this  paper,  we  employ  computational  fluid  dynamic  techniques 
in  order  to  gain  a quantitative  understanding  of  the  physical  trends  in  hybrid  rocket 
combustors. 

The  computational  modeling  is  tailored  to  ongoing  experiments  at  Penn  State  that 
employ  a 2D  slab-burner  configuration.  The  co-ordinated  computational/experimental 
effort  enables  model  validation  while  providing  an  understanding  of  the  experimental 
observations.  Computations  to  date  have  included  the  full-length  geometry  with  and 
without  the  aft-nozzle  section  as  well  as  shorter-length  domains  for  extensive  parametric 
characterization.  HTPB  is  used  as  the  fuel  with  1,3  butadiene  being  taken  as  the 
gaseous  product  of  the  pyrolysis.  Pure  gaseous  oxygen  is  taken  as  the  oxidizer.  The 
fuel  regression  rate  is  specified  using  an  Arrhenius  rate  reaction,  while  the  fuel  surface 
temperature  is  given  by  an  energy  balance  involving  gas-phase  convection  and  radiation 
as  well  as  thermal  conduction  in  the  solid-phase.  For  the  gas-phase  combustion,  a two- 
step  global  reaction  set  is  used.  The  standard  k — e model  is  used  for  turbulence  closure. 
Radiation  is  presently  treated  using  a simple  diffusion  approximation  which  is  valid  for 
large  optical  path  lengths,  representative  of  radiation  from  soot  particles. 

Computational  results  are  obtained  to  determine  the  trends  in  the  fuel  burning 
or  regression  rates  as  a function  of  the  head-end  oxidizer  mass  flux,  G = peUe , and 
the  chamber  pressure.  Furthermore,  computations  of  the  full  slab-burner  configuration 
have  also  been  obtained  for  various  stages  of  the  burn.  Comparisons  with  available 
experimental  data  from  small-scale  tests  conducted  by  General  Dynamics-Thiokol- 
Rocketdyne  suggest  reasonable  agreement  in  the  predicted  regression  rates.  Future  work 
will  include:  (1)  a model  for  soot  generation  in  the  flame  for  more  quantitative  radiative 
transfer  modeling,  (2)  parametric  study  of  combustion  efficiency  and  (3)  transient 
calculations  to  help  determine  the  possible  mechanisms  responsible  for  combustion 
instability  in  hybrid  rocket  motors. 
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Combustion  Processes  in  Hybrid  Rocket  Engines 


S.  Venkateswaran  and  C.  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Presented  at 

The  13th  Workshop  for  CFD  Applications  in 
Rocket  Propulsion,  MSFC,  April  fit5»27,  1995 


Introduction 


© Advantages  of  Hybrid  Propulsion 

— Reduced  Cost 

— Safety 

— Improved  Reliability 

— Thrust  Tailoring 

— Environmentally  Friendly 

• Hybrids  Development 

— Intermittent  Testing  Since  60’ s 

— JIRAD 
— AMROC 

— France  & Japan 

© Small-Scale  Testing 
— JPL/ Strand  et  al. 

— ONERA 
— UAH 
— Penn  State 
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Research  Issues 


• Characterization  of  Fuel  Surface  Regression 

— Fuel  Pyrolysis  and  Surface  Chemistry 

— Heat  Fluxes  - Convection  and  Radiation 

• Combustion  Efficiency 

• Combustion  Instability 

• Modeling  Issue: 

— Boundary  Layer  vs.  Navier-Stokes 
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Expe 


Experimental  Configuration 


• Test  Conditions 

— Fuel  - HTPB 

— Oxidizer  - GOX 

— Pressures  - 300  to  900  psi 

— GOX  Flow  Rates  - 0.2  to  0.8  ibm/s 

©•! 

— GOX  Mass  Flux  (g  = pu)  to  0.5  ibm/in2- 


1320 


Physical  Modeling 


Gas-Phase  Navier- Stokes  Equations 

— Standard  k-e  Model 

Gas-Phase/ Combustion  Model: 

— Butadiene — Product  of  Pyrolysis 

— Two-Step  Global  Kinetics  Model 

C4H6  + 3.5 02  — ► 4 CO  + 3 H20 


CO  + 0.5 02  < — ♦ C02 

Solid-Phase/Pyrolysis: 

— Arrhenius  Pyrolysis  Rate 


Solid/Gas  Coupling 


• Surface  Mass  Balance 

pv  = - psrh 

• Surface  Energy  Balance 

dT  dYi  , f dT  \ , 

~X~dy  + Qrad  + pvh  ~ 2_,  pVim~d^hi  ~ _A<1,  9^  _ P,r‘hs 
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Radiation  Modeling 


• Gaseous  Molecular  Radiation 
— Optically  Thin  Approximation 


4 ckiiTh 


/n  _ M i,j  r- 

Wrad,k  — / y-  j 

■ ■ •-'i  j 


Particulate  (Soot)  Radiation 
— Optically  Thick  Approximation 


Qrad,k  — 


dT 

dy 


5 where  xR  = |tt  £t3 


1323 


Representative  Solution 
Grid  Geometry 


_j T 1 — r- 1 

.00  .30  .60 

Temperature  Contours 
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Representative  Solution 
Axial  Velocity 


Mach  Number  Contours 
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Represent 
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Solutions 


Start  of  Bum 


Midway  Through  Bum 
End  of  Bum 


MASS  FLUX  G 


Representative  Results 
Surface  Regression  Rate 
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Param 

Different 


Start 


.00 


.20 


s 

urn 


3500  K 


Parametric  Studies 
Different  Stages  in  Burn 

With  Radiation/Optically  Thick 
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Parametric  Studies 
Effect  of  GOX  Flow  Rate 


Temperature  Contours 

3500  K 


.00  .20  .40  .60 


Temperature  Profiles 


x=0.3  ro  x=0.45  m x=0.6  ^ 
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Parametric  Characterization  of 
Fuel  Surface  Regression 


• Wffc.  &>&****- 

© Cu^tx4  ^oM2>  W/d 
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Conclusions 


• Navier- Stokes  Analysis  of  Hybrid  Motor 

— Planar  Slab  Burner  Configuration 

— Arrhenius-Rate  for  Pyrolysis 

— Global  Chemistry 

— Turbulence  Model 

— ‘Thick/Thin’  Radiation  Model 

• Computational  Results 

— Parametric  Characterization 

— Fuel  Surface  Temperatures  900  to  1100  K 

— Regression  Rates  of  0.01  to  0.07  in/s 

— Radiative  Fluxes  - Significant  Contribution 

© Ongoing/Future  Work: 

— Radiation  Properties  - Soot  Concentration 

— Combustion  Efficiency  - Downstream  Mixing 

— Combustion  Instabihty  - Transient  Calculations 
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CFD  ANALYSIS  OF  THE  24-INCH  JIRAD  HYBRID  ROCKET  MOTOR 

Pak-Yan  Liang,  Ronald  Ungewitter,  Scott  Claflin 
CFD  Technology  Center 
Rocketdyne  Division,  Rockwell  International 
Mail  Code  IB39,  6633  Canoga  Avenue 
Canoga  Park,  CA  91303 

ABSTRACT 

A series  of  multispecies,  multiphase  CFD  analyses  of  the  24-inch  diameter  joint 
government/industry  IR&D  (JIRAD)  hybrid  rocket  motor  is  described.  The  24-inch 
JIRAD  hybrid  motor  operates  by  injection  of  liquid  oxygen  (LOX)  into  a vaporization 
plenum  chamber  upstream  of  ports  in  the  hydroxyl-terminated  polybutadiene  (HTPB) 
solid  fuel.  Injector  spray  pattern  had  a strong  influence  on  combustion  stability  of  the 
JIRAD  motor  so  a CFD  study  was  initiated  to  define  the  injector-end  flow  field  under 
different  oxidizer  spray  patterns  and  operating  conditions.  By  using  CFD  to  gain  a 
clear  picture  of  the  flow  field  and  temperature  distribution  within  the  JIRAD  motor,  it  is 
hoped  that  the  fundamental  mechanisms  of  hybrid  combustion"  instability  may  be 
identified  and  then  suppressed  by  simple  alterations  to  the  oxidizer  injection 
parameters  such  as  injection  angle  and  velocity. 

The  simulations  in  this  study  were  carried  out  using  the  GALACSY  (General 
ALgorithm  for  Analysis  of  Combustion  SYstems)  multiphase  combustion  codes. 
GALACSY  consists  of  a comprehensive  set  of  droplet  dynamic  submodels 
(atomization,  evaporation,  etc.)  and  a computationally  efficient  hydrocarbon  chemistry 
package  built  around  a robust  Navier-Stokes  solver  optimized  for  low  Mach  number 
flows.  Lagrangian  tracking  of  dispersed  particles  describes  a closely  coupled  spray 
phase. 

The  CFD  cases  described  in  this  paper  represent  various  levels  of  simplification 
of  the  problem.  They  include:  (A)  gaseous  oxygen  with  noncombusting  fuel  vapor 
blowing  off  the  walls  at  various  oxidizer  injection  angles  and  velocities,  (B)  gaseous 
oxygen  with  combusting  fuel  vapor  blowing  off  the  walls,  and  (C)  liquid  oxygen  with 
combusting  fuel  vapor  blowing  off  the  walls.  The  study  used  an  axisymmetric  model 
and  the  results  indicate  that  the  injector  design  significantly  effects  the  flow  field  in  the 
injector-end  of  the  motor.  Markedly  different  recirculation  patterns  are  observed  in  the 
vaporization  chamber  as  oxygen  velocity  and/or  spray  pattern  is  varied.  The  ability  of 
these  recirculation  patterns  to  stabilize  the  diffusion  flame  above  the  surface  of  the 
solid  fuel  gives  a plausible  explanation  for  the  experimentally  determined  combustion 
stability  characteristics  of  the  JIRAD  motor,  and  suggests  how  combustion  stability  can 
be  assured  by  modifications  to  the  injector  design.  Planned  future  activities  to  the 
submodels  which  allow  for  additional  degree  of  realism  will  be  discussed. 
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HYBRID  ROCKET  MOTOR 


. LIQUID  / GAS  OXIDIZER  INJECTED  INTO  VAPORIZATION  CHAMBER 

. SOLID  FUEL  SUBLIMATES  WHICH  THEN  REACTS  WITH 
OXIDIZER 

. INJECTOR-END  RECIRCULATION  PATTERN  DEEMED  CRITICAL  TO 
COMBUSTION  STABILITY  AND  FLAME  HOLDING  MECHANISMS 
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STEP  1 - NON  REACTING  RESULTS 
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CONDITIONS:  TWO  ZONE:  90x40,  100x16 

GOX,  Tin=811°  K,  Uin=132  m/s 

IMAGE  TRUNCATES  COMPUTATIONAL  DOMAIN  AND  PARTICLE  TRACES  DO  NOT  MOVE  BETWEEN 
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